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1.0  Introduction 

Structural  analysis  of  marine  structures  consists  of  the  following  steps: 

a)  description  of  the  environment; 

b)  modeling  of  the  applied  loads; 

c)  load  combination; 

d)  response  analysis,  where  displacements,  nominal  forces  applied  to  each  structural  mem¬ 
ber,  and  stresses  are  calculated; 

e)  fatigue  analysis,  where  damage  inflicted  by  cyclic  loads  is  calculated. 

The  calculated  load  effects  (demand)  are  compared  to  the  strength  of  the  structure 
(capacity),  in  order  to  assess  its  safety. 

In  design,  the  dimensions  of  the  structural  members  are  determined  based  on  the  re¬ 
quirement  that  the  demand  does  not  exceed  the  capacity. 

Uncertainties  are  always  involved  in  all  the  steps  of  structural  analysis  and  also  in 
evaluating  strength.  These  uncertainties  are  due  to  the  random  character  of  the  loading 
environment  and  the  resulting  loads,  or  due  to  inadequate  knowledge  of  physical  phenom¬ 
ena  associated  with  loads. 

Rational  analysis  and  design  of  marine  structures  requires  consideration  of  all  the  un¬ 
certainties  involved  in  predicting  load  effects.  In  probabilistic  methods,  these  uncertainties 
must  be  quantified  in  order  to  assess  structural  safety.  For  example,  if  a  first  order  sec¬ 
ond  moment  (F.O.S.M.)  method  is  used,  we  need  at  least  the  first  two  moments  of  all 
random  variables  associated  with  load  effects  in  order  to  locate  the  design  point  and  eval¬ 
uate  the  safety  index  (Madsen,  86).  Furthermore,  the  determination  of  the  partial  load 
and  resistance  factors,  in  the  safety  equation  of  a  probabilistic  design  code,  also  requires 
quantification  of  all  uncertainties  (ISSC  85). 

The  development  of  probabilistic  analysis  methods  and  design  codes  increased  the  im¬ 
portance  of  quantifying  uncertainties.  Recent  studies  on  offshore  [Faulkner  83,  ISSC  85, 
Das  86,  Guedes  Soares  82],  as  well  as  ship  structures  [Faulkner  81,  ISSC  85,  Kaplan  84, 
Guedes  Soares  84],  investigated  errors  in  evaluating  loads  and  load  effects.  The  results  of 
these  studies  can  be  used  to  assess  the  relative  importance  of  the  various  types  of  uncer¬ 
tainties.  For  example,  one  of  the  conclusions  drawn  from  Guedes  Soaxes  (82)  for  offshore 
structures,  was  that  the  uncertainty  in  the  lifetime  extreme  wave  height  is  the  most  sig¬ 
nificant.  The  error  in  predicting  the  most  severe  sea  condition  over  the  lifetime  is  one  of 
the  major  ingedients  of  the  above  uncertainty. 
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Probabilistic  theory  and  structural  analysis  have  reached  a  state  of  maturity  but  there 
are  still  gaps  in  the  state  of  knowledge  on  quantifying  loads  and  their  effects.  As  part  of  the 
total  effort  associated  with  rational  ship  and  offshore  study  design  based  on  probabilistic 
methods  of  analysis,  a  project  aimed  at  quantifying  the  uncertainties  in  determining  loads 
and  load  effects  in  marine  structures  was  established  by  the  Ship  Structure  Committee. 
The  effort  in  this  project  concentrates  on  assessing  the  quality  of  procedures  for  predicting 
loading  effects  on  marine  structures  by  quantifying  the  errors  involved.  The  following 
issues  are  emphasized  and  discussed  in  the  present  report  summarizing  the  investigation: 

a)  what  is  the  best  way  to  model  uncertainties? 

b)  what  are  the  differences  between  random  (natural)  and  modeling  (subjective)  uncer¬ 
tainties? 

c)  how  do  modeling  uncertainties  affect  extreme  loads? 

d)  what  are  the  most  important  uncertainties? 

1.1  Objective 

The  objective  of  this  study  is  to  quantify  the  error  in  stress  analyses  of  marine  structures, 
thus  providing  necessary  information  to  establish  safety  criteria  in  design. 

The  work  described  herein  is  intended  to  locate  the  sources  of  error  in  all  the  steps  of 
the  load  effect  prediction  procedure  and  to  provide  quantitative  information  on  all  types 
of  error.  The  most  important  types  of  error  are  identified  and  ranked  in  terms  of  then- 
influence  in  design.  Finally,  strategies  axe  recommended  for  reducing  the  most  important 
uncertainties. 

1.2  Report  outline 

The  information  presented  in  this  report  is  organized  as  follows: 

In  Chapter  2,  we  classify  uncertainties  into  two  categories,  random  (natural)  and  mod¬ 
eling,  and  we  study  the  basic  differences  between  these  two  types.  Emphasis  is  given  to 
the  effect  of  modeling  uncertainties  on  extreme  loads.  Furthermore,  we  review  various 
methods  for  modeling  uncertainties. 

The  uncertainties  involved  in  describing  the  loading  environment  for  ships  and  offshore 
platforms  are  considered  in  Chapter  3.  For  ships,  we  studied  the  uncertainty  in  spectral 
shape  variability,  short  crestedness,  directionality  of  weather  systems,  and  visual  obser¬ 
vations.  These  uncertainties  are  quantified  in  terms  of  their  effects  on  the  short  term 
longitudinal  wave  bending  moment.  For  the  case  of  offshore  platforms,  uncertainties  in 
extreme  wave  heights  are  quantified. 
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Chapter  4  deals  with  uncertainties  in  loads  and  load  effects.  For  the  case  of  ships,  most 
of  the  information  is  on  uncertainties  in  short  and  long  term  Stillwater  and  wave  bending 
moments.  Different  ways  for  modeling  uncertainties  are  compared.  It  is  shown  that  we  can 
dramatically  reduce  the  variability,  if  we  use  the  Guedes  Soares  model  for  uncertainty  and 
distinguish  between  different  types  of  ships,  and  between  hogging  and  sagging.  However, 
although  Guedes  Soares’  idea  for  reducing  uncertainties  is  correct,  we  believe  that  it  has 
not  been  properly  implemented  in  Guedes  Soares  (84)  because,  modeling  uncertainties 
have  not  been  correctly  treated  in  his  study. 

For  offshore  platforms,  we  study  uncertainties  in  extreme  global  loads.  Important  fac¬ 
tors,  such  as  current  velocity  and  marine  fouling  are  also  considered. 

Uncertainties  in  load  combination  procedures  are  considered  in  Chapter  5.  The  results 
from  a  Monte  Carlo  simulation  study  on  the  combination  of  wave  and  slamming  induced 
bending  moments  are  presented.  The  objective  of  this  study  is  to  assess  the  error  associated 
with  this  process  and  its  effect  on  extreme  design  loads. 

Chapter  6  deals  with  uncertainties  in  structural  analysis.  The  errors  associated  with 
finite  element  analysis  (F.E.  A.)  of  ships  and  offshore  platforms,  are  quantified.  Information 
is  also  provided  on  the  errors  in  F.E. A.  of  other  types  of  structures  (automotive  and 
aerospace  structures). 

Uncertainties  in  fatigue  analyses  are  studied  in  Chapter  7.  The  study  is  confined  to 
cumulative  damage  based  approaches.  We  examine  the  contribution  of  the  uncertainties, 
which  are  involved  in  all  steps  of  fatigue  analysis,  to  the  overall  uncertainty  in  fatigue 
damage.  This  allows  to  identify  the  most  critical  uncertainties.  Finally,  it  is  shown  that 
the  effect  of  random  uncertainties  on  the  cumulative  damage  is  neglegible  for  both  ships 
and  offshore  structures. 

The  conclusions  from  this  study  are  summarized  in  Chapter  8.  In  this  chapter,  we  also 
identify  the  most  important  uncertainties  involved  in  all  steps  of  stress  analyses  procedures. 
Random  and  modeling  uncertainties  are  compared  in  terms  of  their  effect  on  lifetime 
extreme  loads  and  on  fatigue  damage. 

In  Chapter  9,  we  recommend  a  5  year  research  program.  The  objective  of  this  program  is 
to  expand  the  state  of  knowledge  and  reduce  the  most  critical  uncertainties.  This  program 
has  been  broken  down  into  a  number  of  tasks,  which  have  been  prioritized.  The  tasks  have 
been  prioritized  in  terms  of  the  relative  importance  of  the  associated  uncertainties,  and 
the  expected  benefits,  risk,  and  cost.  A  step  by  step  procedure  is  presented  for  completing 
each  task. 
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2.0  Types  of  Uncertainties 

In  this  chapter  we  define  two  categories  of  uncertainties,  random  and  modeling,  and  ex¬ 
amine  the  differences  between  them.  We  also  review  various  models  for  such  uncertainties. 
Emphasis  is  given  on  the  effect  of  modeling  uncertainties  in  the  extreme  design  loads. 

2.1  Classification 

Uncertainties  in  analysis  of  marine  structures  can  be  categorized  into  natural  (random) 
and  modeling  types.  The  former  axe  due  to  the  statistical  nature  of  the  loading  environment 
and  the  resulting  loads,  and  they  induce  scatter  in  predictions.  The  latter  are  due  to 
the  imperfect  knowledge  on  various  phenomena,  and  idealizations  and  simplifications  in 
analysis  procedures.  These  uncertainties  introduce  both  bias  and  scatter. 

An  example  of  a  natural  uncertainty  is  that  associated  with  the  wave  elevation  at  a  given 
position  in  the  ocean.  An  example  of  a  modeling  uncertainty,  is  the  error  in  calculating 
the  stresses  in  a  structure,  when  the  applied  loads  are  known.  For  this  case,  the  error  is 
only  due  to  the  assumptions  and  simplifications  in  structural  analysis. 

Modeling  uncertainties  are  information  sensitive,  in  the  sense  that  they  can  be  reduced 
as  the  knowledge  of  the  associated  physical  phenomena  expands,  and  the  mathematical 
models  representing  them  become  more  accurate.  This  is  not  the  case  with  random  un¬ 
certainties  which  do  not  decrease  as  we  gather  more  information  on  fundamental  science, 
but  only  as  we  obtain  more  data. 

Both  random  and  modeling  uncertainties  must  be  quantified  and  accounted  for  in  reli¬ 
ability  analysis  and  development  of  probabilistic  design  codes.  In  the  following  paragraph 
we  explain  why  modeling  uncertainties  should  be  considered. 

The  reliability  of  a  structural  system  depends  on  both  load  and  strength  variables. 
Each  variable  can  be  calculated  with  different  degree  of  accuracy.  For  example,  for  most 
of  the  cases,  the  response  of  an  offshore  platform  to  dead  loads  can  be  evaluated  with 
high  accuracy,  given  that  the  loads  are  known,  as  opposed  to  the  case  of  wave  induced 
response  which  cannot  be  predicted  with  the  same  confidence.  Therefore,  when  assessing 
structural  safety  and  making  design  decisions,  we  must  take  into  account  the  differences 
in  the  confidence  levels  associated  with  each  load  and  strength  variable.  For  example,  in 
a  probabilistic  design  code  for  offshore  structures,  the  load  factor  for  wave  loads  is  larger 
than  that  for  dead  loads,  because  the  modeling  uncertainty  associated  with  the  former 
loads  is  larger.  It  should  be  noted,  however,  that  the  load  factor  associated  with  Stillwater 
bending  moments  on  ships  should  be  large  due  to  the  high  coefficient  of  variation  (COV) 
associated  with  these  loads. 
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2.2  Models  for  modeling  uncertainty 

Ang  and  Cornell  (74)  and  Ditlevsen  (82)  proposed  two  different  methods  for  treat¬ 
ing  modeling  uncertainties.  Ang’s  model  is  for  both  load  and  strength  uncertainties. 
Ditlevsen’s  model  was  proposed  for  uncertainties  associated  with  strength  but  it  can  also 
be  applied  to  load  variables. 

In  the  following  we  present  Ang’s  model. 

Let  X  be  the  actual  value  of  some  quantity  of  interest  and  Xq  be  the  corresponding  value 
specified  by  a  design  code.  Then, 


X  =  BiBijX o,  (2.1) 

where  B/  is  the  ratio  of  the  theoretically  predicted  value  for  this  quantity,  Xp,  and  Xq, 
and  Bjj  is  the  ratio  of  X  and  Xp.  Bj  is  a  measure  of  natural  (random)  variability, 
which  is  also  called  type  /  uncertainty,  and  Bn  is  a  measure  of  modeling  uncertainty. 
The  mean  values  of  random  variables  Bj  and  Bjj,  E(Bj)  and  B(B//),  are  the  biases 
corresponding  to  natural  and  modeling  uncertainties,  respectively.  Assuming  that  the 
random  and  modeling  uncertainties  are  statistically  independent,  and  by  using  a  first 
order  second  moment  (F.O.S.M.)  approximation,  which  is  based  on  a  linear  expansion  of 
the  expression  for  B  about  the  mean  value  of  the  random  variables,  we  can  quantify  the 
total  uncertainty  in  A-  as  follows: 

£(B)  =  B(B/)B(B//),  and 

COVB  =  (COV l  +  COVjji  (2'2) 


where  B  =  B/B//. 

COV  stands  for  the  coefficient  of  variation  of  the  quantity  specified  by  the  subscript. 

As  pointed  out  by  Wirsching  (81),  equations  (2.2)  are  valid  for  small  coefficients  of 
variation  (less  than  0.10)  only.  However,  we  adopted  the  above  approximations  because  of 
the  following  reasons: 

a)  They  have  been  employed  in  almost  all  studies  on  uncertainties  in  analyses  on  marine 
structures  (Kaplan  84,  Guedes  Soares  84,  Bea  89,  Olufsen  90,  Moses  85,  86) 

b)  In  most  cases,  we  do  not  know  the  exact  probability  distribution  of  random  and  modeling 
uncertainties.  This  does  not  allow  to  use  more  accurate  methods,  such  as  the  Advanced 
First  Order  Second  Moment  Method  to  combine  random  and  modeling  uncertainties. 

c)  The  estimates  of  the  average  bias  and  COV  of  various  types  of  uncertainty  are  crude. 
Therefore,  it  is  not  reasonable  to  use  more  accurate  (but  also  more  complex)  methods. 
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Random  variables  B[  and  Bn  are  also  assumed  to  be  independent  of  Xq. 

An  example  of  quantifying  modeling  uncertainties  is  illustrated  in  Figure  2.1,  which 
has  been  extracted  from  Bea  (89).  The  quantity  considered  here  is  the  maximum  annual 
wave  height  in  the  northwest  shelf  of  West  Australia.  The  ratio  of  the  measured  over 
the  predicted  maximum  wave  height  is  shown  in  the  horizontal  axis.  The  maximum  wave 
height  is  predicted  using  a  hindcast  method.  The  vertical  axis  represents  the  probability 
that  the  value  of  the  ratio  is  less  than  some  given  number.  Based  on  the  information 
provided  in  Figure  2.1,  the  mean  of  Bn,  which  represents  modeling  uncertainty,  is  1.1 
and  its  coefficient  of  variation  is  0.13.  This  means  that  on  the  average,  hindcast  methods 
predict  a  value  for  the  annual  maximum  wave  height,  which  is  10%  smaller  than  the  actual 
value.  The  coefficient  of  variation  of  0.13  indicates  that  the  ratio  between  the  actual  and 
predicted  wave  height  ranges  between  0.97  and  1.23  with  probability  0.68. 

A  random  variable,  such  as  the  stress  in  a  particular  structural  member,  is  a  function 
of  other  random  variables,  such  as  the  wave  height  and  the  average  wave  period.  Besides 
the  errors  involved  in  calculating  these  variables,  errors  are  also  involved  in  calculating 
the  stress  given  the  values  of  the  latter  random  variables.  Ang  and  Cornell  (74)  presented 
formulas,  which  are  based  on  F.O.S.M.  concepts,  for  quantifying  the  uncertainty  associated 
with  the  above  errors. 

The  Ditlevsen  model  is  applicable  to  reduced  random  variables  (Madsen,  Krenk  and 
Lind,  86),  which  are  independent  gaussian.  We  can  obtain  these  variables  from  the  original 
ones  by  employing  Rosenblatt  transformation  (Rosenblatt  52).  According  to  Ditlevsen 
(82),  model  uncertainty  can  be  accounted  for  by  the  following  equation, 

X'  =  cX'p  +  b  (2.3) 

where  c  is  a  constant,  and  6  is  a  gaussian  random  variable,  which  is  statistically  independent 
of  X' .  The  prime  indicates  reduced  random  variables. 

Ditlevsen,  and  Ang  and  Cornell  models  are  compared  in  Table  2.1. 

Clearly,  Ditlevsen  model  is  more  general.  The  main  difference  between  these  two  mod¬ 
els,  is  that  Ditlevsen  model  accounts  for  the  statistical  correlation  between  the  error  in 
predicting  the  value  of  a  variable,  e  =  X'  —  X'p,  and  the  value  of  the  variable  itself,  while 
Ang’s  model  assumes  that  the  above  random  variables  are  independent.  This  is  demon¬ 
strated  in  Figure  2.2,  which  is  for  the  special  case  that  Xp  and  X  are  lognormal.  The 
value  of  X',  which  is  equal  to  InX,  is  plotted  there  as  a  function  of  Xp.  The  average  of 
X'  or  InX,  as  well  as  regions  corresponding  to  this  average  ±  one  standard  deviation,  are 
plotted  in  this  figure.  It  is  observed  that  the  error  between  actual  values  and  predictions 
for  InX,  which  is  represented  by  the  width  of  the  shaded  region,  is  independent  of  lnXp 
for  the  case  of  Ang’s  model. 
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Note:  X  is  gaussian 

Fig.  2.2  Comparison  Between  Ang  and  Ditlevsen  Models  for 
Modeling  Uncertainty 
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Table  2.1.  Comparison  Between  The  Ang  and  Cornell  and 
Ditlevsen  Models  for  Modeling  Uncertainty 


Characteristic 

Equation: 

Space  in  which  model  is 
applicable: 

Relation  between 
statistics  of  actual 
and  predicted  values: 

Correlation  between  error 
and  predicted  value 


Ditlevsen 

X'  =  cX‘p  +  b 
c:  constant 
b:  random  variable 
independent  of  X'p 

Reduced 

E(X ')  =  cE(Xp)  +  E(b) 
o2x,  =  <?o2x>p  +  a2 

PcX'  =  i  P  . 


Ang 

X  =  BXP 
B:  random  variable 
independent  of  Xp 

Physical 

E(X)  =  E(B)E(XP) 

COVx  =  (COVj  +  COVXp)^ 

P*xp  =  0 


Notes: 


a)  e  denotes  the  error  between  prediction  and  measurements,  i.e.  e  =  X'  -  Xp  (reduced 
space),  or  e  =  X  —  Xp  (physical  space). 


b)  ptx P  denotes  the  correlation  between  e  and  Xp. 


c)  For  the  special  case  of  Xp  and  X  being  lognormal,  Ditlevsen’s  model  reduces  to  Ang’s 
model  for  c  =  1. 

Although  Ang’s  model  is  not  as  general  as  the  Ditlevsen  model,  it  is  preferable,  because  it 
is  simpler.  It  requires  less  information  in  order  to  determine  the  statistics  of  its  parameters, 
and  it  is  very  convenient  to  use  for  the  case  that  the  variables  involved  are  lognormal. 
Moreover,  it  is  expected  that  random  variable  B  is  lognormally  distributed,  for  most  cases, 
because  it  is  usually  the  product  of  several  random  variables.  (Central  Limit  Theorem). 
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2.3  Effect  of  modeling  uncertainties  on  lifetime  extreme  loads 

In  both  ships  and  offshore  platforms,  it  is  important  to  distinguish  between  natural  and 
modeling  uncertainties,  and  the  ways  by  which  they  affect  the  maximum  lifetime  loads  and 
load  effects. 

In  contrast  to  random  uncertainties,  modeling  uncertainties  in  extreme  loads  or  load 
effects  do  not  decrease  with  the  length  of  the  return  period  increasing.  Indeed,  these 
uncertainties  are  systematic.  Consequently,  the  modeling  errors  corresponding  to  two  or 
more  load  applications  are  perfectly  correlated.  Therefore,  the  modeling  error  correspond¬ 
ing  to  the  maximum  of  these  loads  does  not  reduce  with  the  number  of  load  applications 
increasing,  as  it  is  the  case  for  independent  or  weakly  correlated  errors.  Therefore,  uncer¬ 
tainties  in  lifetime  loads  may  be  grossly  underestimated  if  we  treat  modeling  uncertainties 
as  random. 

We  calculate  uncertainties  in  the  extreme  value  of  some  quantity  X  according  to  the 
following  rule:  Let  Xi, ...,  Xn  be  n  independent  samples  from  a  random  variable  and  X^ 
be  their  maximum  value,  i.e. 


X(n)  =max  (Xu...,Xn)  (2.4) 

Then,  the  coefficient  of  variation  of  the  maximum  is: 

COVx(n)  =  -»-  COV/(n) ) *  (2.5) 

where  COV/*(n)  is  the  coefficient  of  variation  corresponding  to  natural  uncertainties  in  the 
maximum  X<"). 

Equation  (2.5)  correctly  implies  that  the  two  types  of  uncertainty,  natural  (random)  and 
modeling,  must  be  treated  differently  when  studying  the  uncertainty  in  the  extreme  value 
of  some  load  or  load  effect.  Furthermore,  the  contribution  of  modeling  uncertainties  to  the 
uncertainty  in  the  maximum  value,  X^n\  does  not  decrease  as  the  number  of  samples,  n, 
increases. 

Olufsen  and  Bea  (90),  and  Bea  (89),  have  concluded  in  their  work  that  uncertainties  in 
maximum  design  loads  and  load  effects  have  been  seriously  underestimated  in  the  recently 
released  API  -  PRAC  22  design  code  for  offshore  platforms.  It  is  remarkable  that  the 
coefficients  of  variation  of  extreme  global  loads,  which  were  derived  from  their  studies,  are 
almost  100%  larger  than  those  used  by  the  developers  of  the  API  code.  This  is  attributed 
to  the  way  in  which  uncertainties  were  treated  in  developing  this  code. 

In  this  report,  we  have  estimated  uncertainties  in  extreme  loads  by  employing  eq.  (2.5) 
for  both  ships  and  offshore  structures.  We  have  demonstrated  that,  for  most  applications, 
this  equation  yields  significantly  larger  uncertainties  than  those  reported  by  some  authors 
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including  the  developers  of  the  API-PRAC  22  design  code.  We  believe  that,  eq.  (2.5)  is 
more  appropriate  than  the  approach  which  is  used  in  the  API-PRAC  22  design  code,  be¬ 
cause  the  former  treats  uncertainties  in  a  more  realistic  way  than  the  latter.  Furthermore, 
we  have  shown  that  the  effect  of  random  uncertainties  on  fatigue  life  is  phased  out  while 
the  effect  of  modeling  uncertainties  does  not  decrease,  with  the  length  of  the  exposure  time 
increasing.  As  a  result,  uncertainties  in  fatigue  life  are  only  due  to  modeling  uncertainties. 
Random  uncertainties  do  not  affect  fatigue  reliability.  The  above  observations  demonstrate 
the  importance  of  modeling  uncertainties. 
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3.0  Loading  Environment 

The  first  step  in  stress  analysis  of  a  marine  structure  is  to  model  the  loading  environ¬ 
ment.  In  this  chapter,  we  study  environmental  uncertainties  associated  with  wave  loads 
on  ships  and  offshore  platforms.  For  offshore  platforms,  we  provide  information  on  other 
types  of  loads  such  as  those  due  to  current  velocity. 

3.1  Ships 

Guedes  Soares  (84)  and  Kaplan  (84)  studied  uncertainties  in  ship  loads  and  load  effects. 
Guedes  Soares  (84)  was  the  first  to  employ  linear  models  relating  the  bias  of  some  quantity 
associated  with  loads  or  load  effects,  to  various  parameters,  such  as  the  significant  wave 
height  and  period.  The  coefficients  in  these  linear  models  were  calculated  by  regression  on 
measured  and  predicted  loads  or  load  effects.  All  other  studies  on  uncertainties  assume  that 
the  bias  is  independent  of  all  parameters.  The  main  advantage  of  Guedes  Soares’  approach 
is  that  it  yields  smaller  coefficients  of  variation  'or  loads  and  load  effects.  Furthermore, 
the  coefficient  of  variation  decreases  when  the  number  of  parameters,  which  are  involved 
in  the  linear  model  for  the  bias,  is  increasing. 

In  this  section  we  examine  uncertainties  in  environmental  description  and  their  effects 
on  short  term  loads  and  load  effects. 

The  following  uncertainties  are  involved  in  modeling  the  loading  environment  for  the 
case  of  ships: 

•  uncertainty  in  wave  spectra, 

•  uncertainty  due  to  short  crestedness, 

•  uncertainty  due  to  directionality  of  the  weather  systems, 

•  uncertainty  in  visual  observations  of  wave  heights, 

•  uncertainty  in  wave  heights  due  to  correlation  of  subsequent  wave  peaks,  and  broad¬ 
bandness  of  the  wave  spectrum, 

•  uncertainties  due  to  heavy  weather  countermeasures. 

Guedes  Soares  (84)  and  Kaplan  (84)  studied  the  uncertainty  in  wave  spectra. 

Guedes  Soares  (84)  separated  the  uncertainty  in  wave  spectra  into  three  categories: 
statistical,  fundamental  and  model  uncertainty.  Statistical  uncertainty  is  associated  with 
errors  in  estimating  a  spectrum  from  experimental  measurements.  Fundamental  uncer- 
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tainty,  which  is  also  called  spectral  shape  variability,  is  due  to  the  variation  of  the  shape  of 
the  spectrum  for  a  specified  sea  state.  Finally,  model  uncertainty  is  due  to  the  discrepancy 
between  the  theoretical  spectrum,  which  is  used  to  describe  a  sea  state,  and  the  average 
of  all  actual  spectra  for  that  sea  state. 

Statistical  uncertainty  was  modeled  by  assuming  that  the  discrepancy  between  the  ordi¬ 
nates  of  the  actual  and  the  theoretical  spectrum  can  be  represented  by  a  gaussian  random 
variable.  The  statistics  of  this  variable  are  the  same  for  all  frequencies,  and  the  random 
variables  representing  the  discrepancies  between  spectral  ordinates  at  different  frequencies 
are  statistically  independent.  This  assumption  was  based  on  Haver  and  Moan’s  (83)  re¬ 
sults.  The  standard  deviation  of  this  variable  for  a  particular  frequency  was  taken  equal 
to  0.525  multiplied  by  the  corresponding  spectral  ordinate  at  that  frequency. 

Guedes  Soares  (84)  estimated  the  bias  and  the  COV  of  the  mean  square  wave  bending 
moment  due  to  statistical  and  fundamental  uncertainties,  for  different  significant  wave 
heights  and  heading  angles.  He  found  that  the  effect  of  statistical  uncertainties  is  negligible 
compared  to  that  of  fundamental  uncertainty  (spectral  shape  variability).  Thus,  the  results 
provided  by  Guedes  Soares  correspond  to  fundamental  uncertainty.  The  effect  of  the 
uncertainty  in  wave  spectra  on  the  mean  square  wave  bending  moment,  for  a  significant 
wave  height  of  10.0  m  and  for  a  Frcude  number  equal  to  0.15,  is  shown  in  Figure  3.1.  This 
figure  has  been  extracted  from  Guedes  Soares  (84).  The  ISSC  spectrum  has  been  used  to 
calculate  the  wave  bending  moment.  Two  cases  of  ships  with  lengths  150m  and  350m  are 
considered.  Results  are  provided  for  different  heading  angles.  The  average  bias  and  the 
COV  over  all  relative  heading  angles  are  also  presented. 


It  is  observed  that  the  error  due  to  uncertainty  in  spectral  shape  is  large  for  small  average 
sea  periods,  for  which  case  the  standard  method  using  the  ISSC  spectrum  is  conservative. 
For  large  average  sea  periods,  (T  >  9sec),  the  error  is  considerably  lower,  with  a  bias 
roughly  equal  to  1.0  and  a  COV  of  0.05. 


An  alternative  simpler  approach,  which  is  called  the  simplified  method,  was  also  used 
by  Guedes  Soares  .  In  the  simplified  method,  it  was  assumed  that  the  bias  depends  only 
on  the  significant  wave  height,  H3.  The  results  of  the  simplified  method  are  summarized 
by  the  following  equations: 


for  ( L  <  250m  or  H„  >  5m) 
0.2 Hs  for  (L  >  250  or  Ha  <  5m) 


(3.1) 


and  COVb.  =  0.1. 


B3  stands  for  the  ratio  between  measured  and  predicted  mean  square  bending  moment, 
E(-)  is  the  mean  value  of  the  quantity  in  parenthesis.  Clearly,  for  long  ships  (L  >  250m) 
and  for  small  wave  heights  ( H„  <  5m),  spectral  shape  variability  leads  to  underestimation 
of  the  wave  bending  moment. 
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Kaplan  (84)  also  studied  the  effect  of  spectral  shape  variability  on  the  root  mean  square 
(rms)  of  the  longitudinal  bending  moment.  He  used  a  set  of  measured  spectra  obtained  by 
Hoffman  and  Miles  (76)  and  evaluated  the  resulting  COV  of  the  rms  bending  moment.  His 
results  show  the  same  trends  as  those  by  Guedes  Soares.  More  specifically,  the  COV  for 
high  sea  states  (Hs  >  10m)  is  significantly  lower  than  that  for  low  significant  wave  heights. 
Moreover,  the  COV  is  lower  for  head  seas  than  for  beam  seas.  Kaplan’s  conclusion  was  that 
the  COV  of  the  rms  bending  moment  due  to  spectral  shape  variability  ranges  between  0. 10 
to  0.20  with  an  average  value  of  0.15.  This  value  is  almost  identical  with  that  estimated 
by  Lewis  (67). 

Uncertainties  are  introduced  by  the  directional  spreading  of  the  wave  spectrum.  Anal¬ 
ysis  procedures,  which  are  recommended  by  all  classification  societies,  use  spreading  func¬ 
tions  to  account  for  directional  spreading.  However,  this  approach  involves  simplifications, 
which  introduce  errors  in  the  calculated  loads.  More  specifically,  it  is  assumed  that  the 
energy  distribution  in  different  directions  is  independent  of  the  significant  wave  height  and 
the  frequency.  This  is  an  unrealistic  assumption,  because  directional  spreading  decreases 
with  the  significant  wave  height  increasing.  Moreover,  the  angular  distribution  of  wave 
energy  is  very  narrow  for  the  frequency  components  near  the  spectral  peak. 

Guedes  Soares  (84)  evaluated  the  bias  and  the  coefficient  of  variation  of  the  mean  square 
bending  moment  due  to  the  effect  of  directional  spreading.  More  specifically,  he  used  the 
following  equation  for  bias, 

=  £  (3.2) 

where  R  is  the  mean  square  bending  moment  calculated  by  using  a  directionality  function, 
which  accounts  for  the  dependency  of  spreading  on  the  significant  wave  height.  Rs  is 
the  mean  square  response  which  is  calculated  by  using  the  directional  spreading  function 
which  is  recommended  by  classification  societies.  The  latter  method  for  calculating  Rs 
will  be  called  standard  method.  The  results  are  plotted  in  Figure  3.2.  Here  the  average 
bias  E(Bac )  is  plotted  against  the  average  sea  period  and  the  heading  angle  for  various 
significant  wave  heights  and  ship  lengths. 

It  is  observed  that  the  average  bias  is  small,  (less  than  0.5),  for  beam  seas,  and  it 
exceeds  1.0  for  head  and  following  seas.  This  means  that  standard  methods  overestimate 
bending  moments  for  beam  seas,  while  they  are  unconservative  for  head  and  following  seas. 
Furthermore,  the  bias  decreases  with  the  significant  wave  height  increasing. 

The  above  observations  can  be  explained  as  follows:  The  wave  energy  spreading,  which 
is  assumed  by  standard  methods,  is  wider  than  the  actual  one.  Moreover,  waves  with 
directions  corresponding  to  heading  angles  of  180°  (head  seas)  induce  considerably  larger 
bending  moments  than  those  corresponding  to  heading  angles  in  the  range  of  90°  (beam 
seas).  Therefore,  standard  methods  underestimate  bending  moments  for  head  seas,  while 
they  overestimate  them  for  beam  seas.  This  can  be  observed  in  Figure  3.2.  The  bias  for 
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Induced  Bending  Moment 


head  seas  increases  for  high  sea  states  because  directional  spreading  decreases  with  the 
significant  wave  height  increasing.  The  opposite  is  true  for  beam  seas.  It  appears  that,  for 
the  cases  studied  by  Guedes  Soares,  the  latter  effect  dominates  and  that  the  average  bias 
decreases  with  significant  wave  height  increasing. 

Finally,  Guedes  Soares  used  a  simplified  approach  for  short  crestedness  bias,  where  the 
bias  was  assumed  to  be  only  a  function  of  the  significant  wave  height.  The  results  of  his 
approach  are  summarized  by  the  following  equations, 

E{Bac)  =  1.  -  0.0077 H,  (3.3) 

COVBte  =  0.05 


Some  refined  formulations  of  the  wave  spreading  function,  which  account  for  the  de¬ 
pendence  of  the  wave  energy  spreading  on  frequency,  have  been  proposed.  However,  it 
is  impractical  to  use  these  formulations  to  calculate  design  loads  because  the  required 
computational  cost  is  too  high.  Moreover,  it  is  very  difficult  to  estimate  some  of  the  pa¬ 
rameters  involved  in  the  spreading  function  models.  Therefore,  we  did  not  estimate  the 
uncertainties  associated  with  these  formulations. 

Due  to  the  systematic  directional  character  of  meteorological  systems,  the  distribution 
of  relative  headings  encountered  by  a  ship  which  always  travels  in  a  certain  route,  is 
nonuniform.  For  example,  a  ship  sailing  eastbound  and  westbound  in  the  Atlantic  Ocean 
is  more  likely  to  encounter  head  than  beam  seas,  as  it  was  found  from  analysis  of  data 
on  ship  headings  (Guedes  Soares  84).  On  the  other  hand,  with  the  exception  of  Bureau 
Veritas,  methods  recommended  by  classification  societies  assume  that  the  distribution 
of  headings  is  uniform  (Liu,  Chen  and  Lee  81).  This  might  introduce  considerable  bias 
because  wave  induced  bending  moments  are  sensitive  to  relative  heading  angles. 

Guedes  Soares  (84)  calculated  the  bias  in  the  mean  square  bending  moment  due  to 
directionality,  for  ships  crossing  North  Atlantic  by  sailing  eastbound  or  westbound.  The 
resulting  bias  is  plotted  in  Figure  3.3  as  a  function  of  the  significant  wave  height  for  various 
wave  periods.  It  is  observed  that  the  standard  method,  which  assumes  that  the  distribution 
of  relative  headings  is  uniform,  is  unconservative.  Moreover,  the  bias  increases  with  the 
significant  wave  height  increasing.  This  can  be  explained  as  follows.  Standard  methods 
assume  that  the  distribution  of  relative  headings  is  uniform,  while,  in  reality,  ships  are 
more  likely  to  encounter  head  seas  than  beam  seas.  Furthermore,  head  seas  induce  higher 
bending  moments  than  beam  seas.  Consequently,  the  actual  bending  moment  is  larger  than 
that  predicted  by  standard  methods.  Moreover,  the  directional  spreading  of  wave  energy 
decreases  with  the  significant  wave  height  increasing,  which  makes  directional  effects  most 
significant.  Thus,  the  bias  due  to  directionality  is  larger  for  high  sea  states. 

The  directional  bias  was  also  calculated  by  Guedes  Soares  by  using  a  linear  model  which 
includes  Ha  as  a  parameter.  For  this  case,  directional  bias,  B<*,  is  given  by  the  following 
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equation, 


(3.3) 


E(Bd)  =  0.981  +  0.018H,. 

The  COV  is  equal  to  0. 10. 

Uncertainty  in  visual  observations  of  wave  heights  also  introduces  uncertainty  in  loads. 
It  has  been  found  that  visual  observations  tend  to  underestimate  small  waves  and  over¬ 
estimate  large  ones.  Guedes  Soares  (84)  recommended  a  linear  model  for  obtaining  the 
significant  wave  height  from  the  visually  observed  wave  height.  This  equation  gives  slightly 
different  results  from  those  obtained  by  the  equation  recommended  by  the  International 
Ship  Structures  Congress  (Hogben  67).  Using  Guedes  Soares  equations,  we  derived  the 
following  equation  for  bias  due  to  uncertaintly  visual  observations, 


E(B„)  = 


0.75  Hs 
Ha  —  2.33  ’ 


(3.4) 


The  COV  was  found  by  Guedes  Soares  equal  to  0.17. 

It  is  important  to  note  that  the  above  values  refer  to  the  error  in  the  significant  wave 
height.  The  corresponding  bias  of  the  mean  square  bending  moment  is  E2(Bu)  and  the 
COV  is  0.34. 


The  effect  of  the  correlation  between  subsequent  wave  peaks  on  the  maximum  lifetime 
wave  height  has  been  studied  by  Naess  (82,83).  Naess  showed  that  the  assumption  of  inde¬ 
pendent  peaks  leads  to  conservative  estimates  of  the  maximum  wave  height.  Dalzell  (89) 
reached  the  same  conclusion  but  he  maintains  that,  for  long  return  periods  corresponding 
to  real  life  applications,  the  error  in  the  expected  maximum  wave  height  is  less  than  10%. 

In  high  sea  states,  the  ship  operator  might  adjust  the  course  in  order  to  reduce  the 
risk  of  capsizing.  For  such  cases,  the  relative  heading  angle  becomes  180°  which  is  the 
direction  corresponding  to  the  largest  wave  bending  moments.  Consequently,  the  effect  of 
countermeasures  in  severe  weather  conditions  is  to  increase  the  wave  bending  moments. 
Guedes  Soares  (90)  performed  a  Monte  Carlo  based  simulation  study  to  assess  the  error 
from  neglecting  the  effect  of  heavy  weather  countermeasures.  He  found  that  the  resulting 
average  bias  is  greater  than  one,  and  that  it  increases  with  the  significant  wave  height 
increasing.  Both  these  results  should  have  been  expected  since  by  neglecting  weather 
countermeasures  we  underestimate  load  effects.  For  a  small  ship  with  length  135m  the 
bias  ranged  between  1.0  and  1.25,  while  for  a  longer  ship  with  length  200m  the  bias  was 
almost  1.0.  The  reason  is  that  heavy  weather  countermeasures  are  very  rarely  taken  for 
large  ships. 

Table  3.1  summarizes  the  uncertainties  involved  in  modeling  the  loading  environment 
for  ships. 
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'Dible  3.1.  Uncertainties  in  Describing  Loading  Environment  for  Ships 


Type  of  uncertainty 


Bias/COV 


Study 


1.  Spectral  shape  variability 


— /0.15 


Kaplan  (84) 


Figure  3.1  (complete  repr.)  Guedes  Soares(84) 

1/0.1  (L  <  250m  or  Ha  >  5m) 

2  -  0.2ffa/0.1(L  >  250m  or  Ha  <  5m) 

(simplified  representation) 


2.  Shortcrestedness 


Figure  3.2  (complete  repr.)  Guedes  Soares(84) 

1  —  0.0077//3/0.05(simplified  repr.) 


3.  Directionality  of  weather 
systems 


Figure  3.3  (complete  repr.)  Guedes  Soares(84) 

0.981  +  0.018i/a/0.10(simplified  repr.) 


4.  Visual  observations 


5.  Heavy  weather  countermeasures 


0.75  Ha 
Ha  -  2.33 


/0.17 


1.0-1.25 


Guedes  Soares(84) 


Guedes  Soares(84), 
Guedes  Soares  (90) 


Notes: 

a)  Kaplan’s  results  refer  to  the  effect  of  spectral  shape  variability  on  rms  wave  induced 
bending  moment.  The  other  results  refer  to  mean  square  bending  moment. 

b)  In  complete  representation,  the  bias  is  considered  as  a  function  of  parameters  such  as  the 
significant  wave  height,  average  wave  period,  and  relative  heading  angle.  In  simplified 
representation,  we  consider  dependence  on  significant  wave  height  only. 
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3.2  Offshore  platforms 

Bea  (89)  investigated  uncertainties  in  global  forces  on  offshore  platforms,  which  are 
located  in  the  northwestern  shelf  of  Western  Australia.  According  to  Bea’s  approach, 
the  maximum  annual  force  is  proportional  to  the  maximum  annual  wave  height  or  its 
square.  The  wave  energy  for  the  particular  platform,  which  was  considered  in  this  study, 
is  out  of  the  range  of  reasonance,  in  which  case  the  period  distribution  is  unimportant. 
Therefore,  the  effects  of  the  loading  environment  were  taken  into  account  by  considering 
the  wave  height  only.  Bea  estimated  both  modeling  and  random  uncertainties  for  the 
expected  annual  maximum  wave  height  by  comparing  hindcast  and  measured  maximum 
wave  heights  in  severe  cyclones.  The  probability  distribution  function  of  the  ratio  of 
measured  and  predicted  annual  maximum  wave  heights,  which  was  used  to  estimate  the 
bias  and  coefficient  of  variation  for  modeling  uncertainties,  is  presented  in  Figure  2.1. 
According  to  this  figure,  the  bias  due  to  modeling  uncertainties  is  lognormally  distributed 
with  median  equal  to  1.1  and  COV  equal  to  0.13. 

The  variability  in  the  expected  severest  sea  condition  over  a  period  of  one  year  is  consid¬ 
ered  in  Figure  3.4,  which  is  extracted  from  Bea  (89).  The  probability  distribution  function 
of  the  annual  maximum  significant  wave  height  is  plotted  in  this  figure.  The  maximum 
wave  height  is  also  considered  in  the  same  figure.  According  to  these  results,  the  COV 
of  the  random  bias  of  the  annual  maximum  wave  height  is  0.30.  Assuming  that  modeling 
and  random  uncertainties  are  independent,  the  COV  due  to  both  types  of  uncertainties  is 
0.33. 

Uncertainty  due  to  the  facts  that,  a)  the  wave  peaks  do  not  follow  the  Rayleigh  distri¬ 
bution,  and  b)  subsequent  wave  peaks  are  correlated  has  been  studied  by  Forristall  (78). 
Based  on  analysis  of  data  from  the  Gulf  of  Mexico,  Forristall  found  that  the  assumption 
associated  with  the  Rayleigh  distribution  and  the  independence  of  wave  peaks  can  over¬ 
estimate  the  wave  heights  by  10  to  20  percent.  Bea  (89)  used  Forristall’s  results,  and 
recommended  an  average  bias  for  the  wave  height  equal  to  0.93  and  a  COV  of  0.08.  Con¬ 
sequently,  the  bias  in  the  mean  square  bending  moment  due  to  this  effect  is  0.86  and  the 
COV  is  0.16. 

Wirsching  (81)  also  studied  the  uncertainty  on  load  and  load  effects  due  to  uncertainties 
in  modeling  the  loading  environment.  His  study  was  confined  to  modeling  uncertainties 
associated  with  wave  loads.  Data  on  uncertainties  were  also  reported  in  Wirsching  (88). 
According  to  Wirsching,  the  average  bias  ranges  between  0.63  to  1.20  and  the  COV  ranges 
between  0.40  to  0.60.  For  offshore  platforms  with  low  natural  frequencies,  for  which  dy¬ 
namic  effects  are  important,  the  largest  component  of  this  uncertainty  is  due  to  variability 
in  the  dominant  period  of  the  spectrum  describing  a  sea  state.  Although  the  variability 
in  the  dominant  period  might  be  moderate  (the  COV  ranges  between  0.10  and  0.20),  its 
effect  on  the  resulting  loads  and  load  effects  is  high.  The  high  sensitivity  of  a  platform’s 
response  to  uncertainties  in  dominant  period  has  also  been  discussed  by  Moses  (85).  Un- 
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certainties  in  the  relative  frequency  of  occurrence  of  each  sea  state,  as  well  as  those  due 
to  directionality  of  weather  systems  have  also  been  accounted  for  by  Wirsching  in  quan¬ 
tifying  uncertainties  in  environmental  description.  The  median  bias  of  the  former  type  of 
uncertainties  ranges  between  0.9  and  1.00  and  the  COV  between  0.1  and  0.3. 

Guedes  Soares  and  Moan  (82)  investigated  uncertainties  in  global  loads  applied  to  a 
cylindrical  pile  in  the  North  Sea.  The  COV’s  for  the  model,  random  and  measurement 
uncertainties  in  the  maximum  wave  height  were  found  equal  to  0.11,  0.08  and  0.09,  respec¬ 
tively.  The  maximum  wave  height  corresponds  to  a  period  of  100  years. 

A  linear  model  was  postulated  for  the  period  corresponding  to  the  maximum  wave 
height,  and  the  model  parameters  were  estimated  by  regression.  The  COV’s  for  the  model, 
random  and  measurement  uncertainty  in  the  wave  period  were  estimated  equal  to  0.10, 
0.05  and  0.10,  respectively.  The  total  COV  for  the  period  is  0.14,  and  the  correlation 
coefficient  between  the  wave  height  and  the  period  is  0.50,  which  implies  a  strong  positive 
correlation  between  these  quantities. 

Guedes  Soares  and  Moan  also  considered  current  loads.  The  COV  of  the  current  velocity 
was  estimated  0.35.  No  information  was  provided  on  the  probability  distribution  of  the 
aforementioned  variables.  Furthermore,  the  error  due  to  shortcrestedness  of  waves  was  not 
examined. 

A  survey  was  conducted  under  the  sponsorship  of  ABS  and  Conoco  in  order  to  quantify 
the  modeling  error  associated  with  design  loads  for  cylinders  and  pontoons  for  tension  leg 
platforms  (Wirsching  82).  The  bias  and  the  COV  for  the  cumulative  effect  of  wind,  current 
and  tide  on  the  platform  loads  were  reported  equal  to  0.95  and  0.12,  respectively.  These 
estimates  were  based  on  the  responses  of  experts  to  a  questionnaire. 
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Table  3.2  summarizes  the  uncertainties  involved  in  modeling  the  loading  environment 
for  offshore  platforms. 

Table  3.2.  Uncertainties  in  Describing  Environment  for  Offshore  Platforms 
Type  of  uncertainty  Bias/COV  Source 


Maximum  wave  height 

1.0/0.30  (Random)1 
1.1/0.13  (Modeling) 

Bea  (89) 

Maximum  wave  height 

-/0.162 3 

Guedes  Soares 
and  Moan  (82) 

Average  period  corresponding  — /0.142,3 

to  maximum  wave  height 

Guedes  Soares 
and  Moan  (82) 

Effect  of  environmental 
uncertainties  on  rms 
stress  in  structural  elements 
of  a  platform 

0.6  -  1. 2/0.4  -  0.64 

Wirsching(81) 

Effect  of  modeling 
uncertainties  in  wind, 
current  and  tide  on  loads 
on  tension  leg  platform 

0.95/0.12 

Wirsching(82) 

Notes: 

1  Results  refer  to  annual  maximum  wave  height.  Bias  is  lognormally  distributed.  Un¬ 
certainties  due  to  the  correlation  of  subsequent  peaks  and  the  Rayleigh  distribution 
assumption  are  included  in  these  results. 

2  Results  refer  to  100  year  maximum  wave  height.  The  error  due  to  short  crestedness  was 
not  examined. 

3  Correlation  coefficient  between  wave  height  and  period  is  0.5. 

4  Bias  is  lognormally  distributed.  The  relation  of  the  period  to  the  response  depends 
upon  the  natural  period  characteristics  of  the  structure.  Uncertainty  in  average  wave 
period  dominates  for  offshore  platforms  with  low  natural  frequencies,  for  which  there  is 
a  substantial  amount  of  wave  energy  in  the  vicinity  of  their  natural  frequency 
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4.0  Loads 

In  this  chapter,  we  study  uncertainties  in  loads  and  load  effects.  For  ships,  we  exam¬ 
ine  loads  applied  to  the  main  girder  as  well  as  hydrodynamic  pressure.  Uncertainties  in 
both  short  and  long  term  predictions  are  quantified.  For  offshore  platforms,  we  quantify 
uncertainties  in  base  shear  and  overturning  moment. 

4.1  Ships 

Here,  we  consider  uncertainties  in  the  most  important  loads  applied  to  ships. 

4.1.1  Stillwater  bending  moments  and  shear  forces 

Guedes  Soares  and  Moan  (88)  statistically  analyzed  Stillwater  bending  moments  and 
shear  forces  for  various  ship  types.  In  their  study,  Stillwater  load  effects  were  assumed  to 
vary  from  voyage  to  voyage  for  a  particular  ship,  from  one  ship  to  another  in  a  particular 
class  of  ships,  and  from  one  class  of  ships  to  another.  The  above  sources  of  variability  can 
be  modeled  as  follows, 


mijk  =  m0  +  mic  +  rrij  +  ti 


(4.1) 


where, 

rriijk  is  the  bending  moment  or  shear  force,  at  the  ith  voyage,  which  is  applied  to  the  jth 
ship,  which  belongs  to  the  kth  class, 


m0  is  the  average  load  effect  for  all  ships, 

m0  +  rrifc  is  the  average  load  effect  of  all  ships  in  the  kth  class, 

m0  +  rrik  +  m j  is  the  average  load  effect  for  the  jth  ship  of  the  kth  class, 


and  €j  represents  the  variation  of  the  load  effect  from  voyage  to  voyage.  Accordingly,  the 
following  variances  can  be  defined, 

a)  variance  of  the  load  effect  for  a  particular  ship:  a\ , 

b)  variance  of  the  load  effect  for  all  ships  in  a  particular  class,  which  is  specified  by  k: 

c)  variance  of  the  load  effect  for  all  ships:  ( a 2  +  cr2  +  o*)1/2. 

The  generality  of  description  increases  from  a)  to  c)  by  accounting  for  all  ships  in  a 
class,  or  by  accounting  for  all  ships  in  all  classes.  Clearly,  the  variance  increases  with 
the  generality  of  description  increasing. 
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Tables  4.1  and  4.2  summarize  the  results  from  statistical  analysis  of  data  on  Stillwater 
bending  moments  and  shear  forces  respectively,  for  seven  types  of  ships.  The  values  in 
these  tables  have  been  normalized  by  the  corresponding  values  which  are  prescribed  by 
classification  societies.  The  average  Stillwater  load  effect,  and  the  variance  of  this  load 
effect  for  one  ship,  and  also  for  all  ships  in  a  given  class,  are  presented  in  Table  4.1.  The 
results  are  based  on  the  analysis  performed  by  Guedes  Soares  and  Moan  (88),  and  all  the 
numbers  are  normalized  by  the  design  values  prescribed  by  classification  societies.  The 
data  used  in  this  analysis  can  be  found  in  Guedes  Soares  and  Moan  (88)  and  in  Guedes 
Soares  (84). 

Table  4.1  Variability  in  Stillwater  bending  moments 


Type  of  ship 

Average  load  effect 
for  all  ships 
within  a  class 

(ra0  4-  mk) 

Standard  deviation 
for  a  particular 
ship 

M 

Standard  deviation 
for  all  ships 
within  a  class 

(al  +  a])^ 

Cargo 

0.50 

0.28 

0.30 

Containership 

0.72 

0.16 

0.20 

Bulk  Carrier 

-0.008 

0.30 

0.38 

OBO 

0.80 

0.30 

0.41 

Chemical  Carrier 

-0.005 

0.22 

0.36 

Ore/Oil  Carrier 

-0.44 

0.22 

0.37 

Tanker 

-0.12 

0.21 

0.44 

Note: 

a)  Positive  bending  moments  correspond  to  hogging,  and  negative  ones  to  sagging. 

b)  the  bending  moments  have  been  normalized  by  dividing  by  the  values  which  are  pre¬ 
scribed  by  classification  societies. 
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Table  4.2  Variability  in  stiUwater  shear  forces 


Type  of  ship 

Average  load  effect 
for  all  ships 
within  a  class 

(m0  +  mk) 

Standard  deviation 
for  a  particular 
ship 

M 

Standard  deviation 
for  all  ships 
within  a  class 

(*e2  +  *?)1/2 

Cargo 

-0.17 

0.20 

0.31 

Containership 

-0.32 

0.10 

0.29 

Bulk  Carrier 

0.04 

0.34 

0.52 

OBO 

0.03 

0.25 

0.49 

Chemical  Carrier 

-0.06 

0.20 

0.30 

Ore/Oil  Carrier 

0.12 

0.25 

0.75 

Tanker 

0.07 

0.20 

0.55 

It  is  observed  that  cargo  and  container  ships  experience  large  hogging  moments.  Tankers 
and  Ore/Oil  Carriers  axe  subjected  to  sagging  moments.  Although  the  average  Stillwater 
bending  moment  is  small  for  tankers,  there  is  a  large  variability  in  this  moment.  This  is 
attributed  to  the  large  variability  of  the  Stillwater  bending  moment  from  one  tanker  to 
another. 

Kaplan  (84)  reported  some  results  on  Stillwater  bending  moments  obtained  from  Akita 
(82).  According  to  his  study,  the  COV  for  containerships  is  0.29,  and  for  tankers  it  is  0.99 
for  ballast,  and  0.52  for  full  load  conditions.  These  values  reflect  variabilities  from  voyage 
to  voyage  and  from  one  ship  to  another  within  a  particular  class.  They  indicate  the  same 
trend  with  Guedes  Soares  results.  Indeed,  the  variability  is  considerably  larger  for  tankers 
than  for  containerships.  However,  Guedes  Soares  reported  a  significantly  larger  variability 
for  tankers  (COV  '-±  3.7 )  compared  to  that  reported  by  Kaplan.  This  discrepancy  might  be 
due  to  the  large  spreading  of  sizes  of  the  tankers  which  were  considered  by  Guedes  Soares. 
The  coefficients  of  variation  reported  for  containerships  are  almost  identical 

The  above  results  are  for  Stillwater  load  effects,  which  occur  at  any  time  instant.  The 
extreme  values  of  these  load  effects  are  also  important.  In  the  following,  we  calculate 
the  first  two  moments  of  the  lifetime  maximum  Stillwater  bending  moments,  by  using  the 
results  from  Tables  4.1  and  4.2. 

The  value  of  a  load  effect  at  a  particular  voyage  is  given  by  eq.  (4.1).  c7  is  assumed  to 
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be  a  zero  mean  gaussian  random  variable  with  standard  deviation  ae.  Furthermore,  the 
values  of  c  for  different  voyages  are  assumed  to  be  independent. 

The  lifetime  extreme  load  effect  for  a  particular  ship  is; 


maxfm^A:}  =  m0  +  mk  +  mj  +  max  ej  (4.2) 

Subscript  i  below  max  indicates  that  the  maximum  value  of  the  load  effect  refers  to  all 
voyages. 

Therefore,  the  standard  deviation  of  the  lifetime  maximum  load  effect  for  all  ships  in  a 
particular  class  is, 

Om  =  {a]  +  ^LXJ1/2  (4-3) 

where  crmax«  denotes  the  standard  deviation  of  the  maximum  load  effect,  maxije,}.  The 
latter  random  variable  follows  the  type  I  asymptotic  extreme  value  probability  dist  ribution 
(Ochi  90), 

-fmaxe(e)  =  exp  (exp  {-aN  (e  -  bN )))  (4.4) 

where, 

N  is  the  number  of  voyages  over  the  lifetime  of  the  ship, 
aN  =  N  f  ( bN ), 


/  (•)  is  the  probability  density  function  followed  by  6j, 


bN  is  the  most  probable  value  of  max  e,,  which  is  equal  to  4*_1 
and  $  (•)  is  the  standard  gaussian  probability  distribution. 

The  mean  and  variance  of  max  e*  are, 


E  (maxej)  =  bff -\ - ,and 

*  on 


_2  = 
max£  6  al 


respectively.  7  is  the  Euler  constant. 


Table  4.3  presents  the  results  for  lifetime  extreme  Stillwater  bending  moments  for  N  — 
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Table  4.3  Uncertainties  in  lifetime  extreme  Stillwater  bending  moments 


Type 
of  ship 

Most  probable 
extreme  bend, 
moment 

Stdr.  dev. 
for  variability 
from  ship  to 
ship  within  a 
particular  class 

(<u) 

Stdr.  dev. 
of  extreme 
bend,  moment 
for  a  par¬ 
ticular  ship 

(^max  e) 

Stdr.  dev. 
of  extreme 
bend,  moment 
for  all  ships 
within  a  class 

(<*.,«+<’]) 1/2 

Cargo 

1.27 

0.11 

0.12 

0.16 

Containership 

1.16 

0.12 

0.07 

0.14 

Bulk  Carrier 

-0.84 

0.24 

0.13 

0.27 

OBO 

1.13 

0.28 

0.13 

0.31 

Chemical  Carrier 

-0.66 

0.29 

0.10 

0.31 

Ore/Oil  Carrier 

-1.04 

0.30 

0.10 

0.32 

Tanker 

-0.70 

0.39 

0.09 

0.40 

It  is  observed  that,  for  the  case  of  Chemical  Carriers,  Ore  Carriers  and  Tankers,  the 
variability  in  lifetime  extreme  Stillwater  bending  moments  is  almost  entirely  due  to  the 
variation  of  bending  moments  from  one  ship  to  another  within  a  particular  class.  Indeed, 
for  this  type  of  ships,  the  variability  from  one  voyage  to  another  is  too  small  to  affect  the 
extreme  bending  moments. 

4.1.2  Uncertainties  in  short  term  vertical  wave  bending  due  to  errors  in  re¬ 
sponse  amplitude  operators 

Kaplan  (84)  compared  model  data  against  theoretical  predictions  of  response  ampli¬ 
tude  operators  for  two  Series  60  ships  (0.70  and  0.80  block  coefficients),  and  also  for  the 
WOLVERINE  STATE.  The  data,  which  can  be  found  in  Kaplan  and  Raff  (72),  cover 
different  speeds  and  headings  in  regular  waves.  Kaplan  calculated  the  rms  of  the  wave 
bending  moment  by  using,  a)  theoretically  calculated  response  amplitude  operators  ob¬ 
tained  from  the  SCORES  seakeeping  computer  code,  and  b)  measured  response  amplitude 
operators.  A  reference  wave  spectrum  was  used,  for  which  the  value  of  the  power  spectral 
density  function  is  constant  with  frequency.  The  bias  due  to  errors  in  response  amplitude 
operators  was  calculated  by  comparing  the  rms  values  of  the  wave  bending  moment,  which 
were  calculated  by  using  experimental  and  theoretical  response  amplitude  operators. 
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Based  on  the  above  approach,  Kaplan  found  that  the  COV  of  the  rms  wave  bending 
moment  is  0.10.  No  information  was  provided  on  the  probability  distribution  of  the  bias 
or  its  average  value. 

Guedes  Soares  (84)  separated  uncertainties  in  response  amplitude  operator  into  those 
due  to  nonlinearities  and  those  due  to  all  the  other  simplifications  and  idealizations.  Ac¬ 
cording  to  his  approach  the  bias  in  the  response  amplitude  operator  is  given  by  the  following 
equation; 

H(u)  =  Bl  BS/h  Hp  (^)  for  any  u >  (4.6) 

where,  Bl,  is  the  bias  due  to  all  uncertainties  except  nonlinearities,  B$  expresses  the  un¬ 
certainty  in  sagging,  Bu  expresses  the  uncertainty  in  hogging,  H  is  the  actual  response 
amplitude  operator,  and  Hp  is  the  value  of  the  response  amplitude  operator  as  it  is  pre¬ 
dicted  by  a  linear  strip  theory  based  method.  Errors  due  to  the  flexibility  of  the  ship  hull 
were  found  to  be  unimportant  except  for  very  long  ( L  >  350m),  fast  ships.  Therefore,  this 
source  of  uncertainty  was  neglected.  Guedes  Soares  examined  the  error  associated  with 
the  Salvensen,  Tuck  and  Faltinsen  (70)  (S.T.F.)  method.  Linear  models  were  postulated 
for  both  Bl,  and  BS/H,  and  the  coefficients  were  found  by  regressing  on  data  from  model 
experiments.  Bl,  was  assumed  to  be  a  function  of  the  relative  heading  angle  a,  the  Froude 
number  V,  and  the  block  coefficient  CB •  The  following  relations  were  found  for  the  bias, 
on  the  basis  of  regression  fits: 

Bl  =  0.00631a  +  1.22V  +  0.657 CB  +  0.064  for  0  <  a  <  90°,  and  (4.7) 

Bl  =  -0.00495a  +  0.42V  +  0.701CB  +  1.28  for  90°  <  a  <  180°. 

The  COV  was  found  equal  to  0.38  for  both  cases. 

An  alternative  approach  was  also  followed,  in  which  the  linear  bias,  Bl,,  was  assumed 
to  be  a  function  of  the  significant  wave  height,  Hs,  only.  For  this  case,  the  bias  is  given 
by  the  following  equation, 

Bl  =  1.22  -  0.005i/s.  (4.8) 

The  COV  was  found  equal  to  0.35.  The  bias  in  eq.  (4.8)  is  defined  as  the  ratio  of 
the  average  values  of  the  measured  and  predicted  response  amplitude  operators  over  all 
heading  angles  and  average  wave  periods. 

The  effect  of  nonlinearities  was  modeled  by  employing  a  linear  model  which  involved 
the  block  coefficient  CB  as  a  parameter.  The  resulting  equations,  which  were  also  derived 
by  regression,  are, 

Bs  —  1.74  —  0.93CB  for  sagging,  and  (4.9) 

Bh  —  0.26  +  0.93CB  for  hogging. 

The  COV  was  found  equal  to  0.12  for  both  equations. 
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The  following  conclusions  can  be  extracted  from  eq’s  (4.7)-(4.9): 

•  S.T.F.  method  is  unconservative,  when  it  is  used  to  predict  sagging  bending  moments. 

•  The  error  of  the  S.T.F.  method  is  larger  for  beam  seas  than  it  is  for  head  and  following 
seas.  For  example,  the  bias,  Bl,  for  a  =  90°,  V  =  0.2,  and  Cb  =  0.8,  is  1.48,  while  it 
is  only  1.03  for  a  =  180°  and  same  V  and  Cb • 

•  S.T.F.  method  underestimates  sagging  and  overestimates  hogging  because  the  linear 
model  does  not  distinguish  between  them. 

•  The  error  of  the  S.T.F.  method  due  to  nonlinearities  is  smaller  for  ships  with  large 
block  coefficients.  This  is  true  because  the  assumption  of  vertical  hull  walls  is  realistic 
for  ships  with  large  Cb- 

•  The  bias,  Bl,  decreases  with  Hs  increasing. 

Although  a  large  portion  of  the  experimental  data  used  by  Kaplan  and  Guedes  Soares 
are  identical,  a  significant  discrepancy  is  observed  between  their  COV’s.  In  our  opinion,  the 
above  discrepancy  should  be  attributed  to  the  way  by  which  uncertainties  were  quantified 
by  Guedes  Soares  (84).  More  specifically,  Guedes  Soares  regressed  on  data  for  the  ratio  of 
measured  and  predicted  response  amplitude  operators  for  various  frequencies.  A  typical  set 
of  data  is  shown  in  Figure  4.1.  Clearly,  this  approach  overestimates  modeling  error,  because 
it  uses  data  from  test  measurements  which  have  been  contaminated  with  experimental 
errors  as  well  as  concentrating  on  individual  frequencies.  A  better  way  to  proceed  is  the 
following, 

a)  postulate  a  linear  model  for  the  rms  bending  moment, 

b)  transform  the  data  on  transfer  function  into  data  on  the  rms  bending  moment  by  using 
some  sea  spectrum  (for  example,  the  ISSC  spectrum)  and  by  integrating  over  frequency, 

c)  regress  on  the  data  from  b),  or  simply  estimate  the  COV  of  the  ratio  of  measured  over 
predicted  rms  bending  moments. 

This  procedure,  which  has  been  followed  by  Kaplan  (84),  removes  the  experimental 
error  as  well  as  the  individual  frequency  sensitivity  by  integrating  over  frequency  in  step 
b).  Therefore,  in  our  opinion,  the  results  obtained  from  this  approach  should  be  more 
realistic. 

4.1.3  Long  term  induced  bending  moments 

Kaplan  (84)  found  that  the  COV  of  the  extreme  lifetime  vertical  bending  moment  is 
0.19.  The  COV  of  random  uncertainties  was  found  0.065.  No  information  was  provided 
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Figure  4.1  Calculated  and  Measured  Response  Amplitude 
Operators  for  a  Containership  (Fn=0.245) 


on  the  probability  distribution  of  the  average  value  of  the  bias.  The  relative  contribution 
of  the  uncertainties  examined  by  Kaplan  is  presented  in  Table  4.4. 

Guedes  Soares  (84)  estimated  uncertainties  in  the  most  probable  extreme  long  term 
vertical  bending  moment  for  the  following  cases, 

a)  tankers  ( Cb  =  0.8) 

b)  tankers  in  sagging, 

c)  containerships  (Cb  =  0.6)  in  hogging, 

d)  containerships  (Cb  =  0.6)  in  sagging, 

e)  hogging  in  any  type  of  ship, 

f)  sagging  in  any  type  os  ship, 

g)  the  type  of  ship  and  bending  moment  (hogging  or  sagging)  is  unknown. 

An  extension  of  Ochi’s  approach  was  used  to  caluclate  the  long  term  extreme  long  term 
vertical  bending  moment  (Ochi  78).  The  results  from  Guedes  Soares  study  are  shown  in 
Table  4.5. 

Clearly,  the  variability  in  load  effects  is  smaller  for  cases  that  the  type  of  ship  and/or  the 
type  of  moment  are  specified  in  the  formulation.  For  example,  the  modeling  bias  for  a  ship 
with  block  coefficient  of  0.8  is  1.13  and  the  COV  is  only  0.04.  If  the  block  coefficient  is  not 
specified,  the  bias  is  1.10  and  the  COV  is  0.15.  This  indicates  that  a  design  code,  which 
distinguishes  between  various  ship  types  and  hull  characteristics  and  specifies  different  load 
and  strength  factors  for  each  case,  allows  the  design  of  more  efficient  ships. 

Another  conclusion  from  Guedes  Soares  study  is  that  theoretical  predictions  are  almost 
always  unconservative.  This  is  primarily  due  to  the  unconservative  errors  of  linear  strip 
theory  in  response  amplitude  operators.  In  particular,  the  error  in  the  sagging  bending 
moments  is  very  large  for  ships  with  small  block  coefficients.  For  example  the  bias  is  1.28 
for  containerships.  The  reason  is  that  nonlinearity  in  response  is  significant  for  these  ships, 
due  to  their  nonvertical  sides.  This  unconservative  error  must  be  accounted  for  in  design 
because  sagging  can  cause  buckling  of  the  deck  plates,  which  is  an  important  failure  mode 
in  ship  hulls. 
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Table  4.4  Relative  Contribution  of  Various  Types  of 
Uncertainties  to  Total  Uncertainties  in  Extreme 
Bending  Moment  (source:  Kaplan  (84)) 


Type  of  Uncertainty 

Contribution 

(%) 

Spectral  shape  variability 

61 

Uncertainty  in  transfer  function 

27 

Random  uncertainty 

12 

Note:  Contribution  is  defined  as  the  square  of  the  ratio  of  the  particular  uncertainty  over 
the  total  uncertainty. 


Table  4.5  Uncertainties  in  Long  Term  Vertical  Wave 
Bending  Moment  (source:  Guedes  Soares  (84)) 


Case 

Modeling 

Uncertainty 

Exp.  Bias/COV 

Tankers  in  hogging 

1.13/0.04 

Tankers  in  sagging 

1.13/0.04 

Containerships  in 
hogging 

0.88/0.05 

Containerships  in 
sagging 

1.28/0.04 

Any  ship  in  hogging 

1.0/0.15 

Any  ship  in  sagging 

1.2/0.08 

Any  ship/hogging 
or  sagging 

1.1/0.15 

Random  Uncertainty  Total  Uncertainty 
COV  Exp.  Bias/COV 


0.07 

1.13/0.08 

0.07 

1.13/0.08 

0.07 

0.88/0.09 

0.07 

1.28/0.08 

0.07 

1.0/0.17 

0.07 

1.2/0.11 

0.07 

1.1/0.17 

27 


Although  Guedes  Soares  and  Kaplan’s  results  are  in  good  agreement,  we  believe  that 
modeling  uncertainties  were  not  treated  properly  by  the  former.  More  specifically,  Guedes 
Soares  assumed  that  modeling  errors  in  mean  square  bending  moments  are  independent 
from  one  sea  state  to  another  or  from  one  heading  to  another  (eq.  4.57,  p.  278  of  Guedes 
Soares  (84)).  This  assumption  is  not  realistic,  because,  as  we  mentioned  in  chapter  2, 
modeling  uncertainties  are  systematic  and  as  such,  they  are  highly  correlated  from  one  sea 
state  to  another  or  from  one  heading  angle  to  another.  Therefore,  the  COV’s  reported  in 
Guedes  Soares  (84)  might  be  lower  than  the  actual  values. 

Faulkner  (81)  reported  the  following  COV’s  for  lifetime  extreme  vertical  bending  mo¬ 
ments: 

a)  Modeling  uncertainties: 

0.15  for  warships 

0.10  for  commercial  ships 

b)  Random  uncertainties: 

0.12  for  both  warships  and  commercial  ships. 

Faulkner  considered  a  SL-7  containership  and  a  large  tanker  in  his  study.  He  found 
that  the  uncertainty  is  larger  for  the  containership  than  for  the  tanker,  which  agrees  with 
Guedes  Soares’  conclusions. 

Finally,  uncertainties  in  both  vertical  and  horizontal  bending  moments  were  considered 
in  ISSC(85),  for  a  tanker  with  length  equal  to  160m.  The  bias  of  both  bending  moments 
was  assumed  to  be  normal  with  a  mean  of  0.95  and  a  COV  of  0.1  for  the  vertical  bending 
moment.  The  bias  and  COV  for  the  horizontal  bending  moment  axe  0.85  and  0.15,  re¬ 
spectively.  The  correlation  coefficient  between  the  two  bending  moments  was  assumed  to 
be  0.70.  Unfortunately,  no  information  was  provided  in  ISSC(85)  on  how  these  numbers 
were  derived.  Moreover,  it  is  mentioned  in  this  report  that  these  numbers  are  simply  crude 
approximations. 

i'he  results  from  the  studies  considered  in  this  section  axe  summarized  in  Table  4.6. 
4.1.4.  Uncertainties  in  hydrodynamic  pressure 

Chen  et  al  compared  theoretically  predicted  hydrodynamic  pressures  on  a  ship  hull 
against  model  tests  results  and  full  scale  measurements.  A  linear  strip  theory  based  com¬ 
puter  code  (ABS/SHIPMOTION)  was  used  to  calculate  pressures.  Measurements  were 
obtained  for  an  SL-7  containership  and  a  Great  Lakes  bulk  carrier.  The  total  hydrody¬ 
namic  pressure,  the  pressure  component  due  to  the  incident  and  diffracted  waves,  and  the 
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pressure  component  arising  from  ship  motions  were  considered  in  this  study.  Model  tests 
were  performed  for  head  seas  at  Froude  numbers  0.15,  0.23  and  0.32  over  a  range  of  ship 
length/wave  length  ratios  from  0.65  to  1.65. 


The  following  are  the  main  conclusions  from  Chen’s  study: 

•  The  calculated  pressures  due  to  ship  motions  correlated  well  with  test  measurements. 

•  Good  agreement  was  also  found  between  predictions  and  measurements  for  the  pressure 
due  to  incident  and  diffracted  waves. 

•  The  agreement  between  predictions  and  measurements  for  the  total  hydrodynamic  pres¬ 
sure  was  good  except  for  the  bow  and  stern  regions.  This  should  have  been  expected 
because  three  dimensional  effects  and  nonlinearities  are  stronger  in  these  regions. 


Table  4.6  Summary  of  Results  On  Uncertainties 
in  Long  Term  Extreme  Bending  Moments. 

Bias/COV 


Source 

Quantity 

Kaplan  (84) 

Guedes  Soares 
(84) 

ISSC  (85) 

Faulkner  (81) 

V.  Bend.  Mom. 
(*i) 

— /0.19 

0.88  -  1.28/ 
0.08-0.17 

0.95/0.1 

— /0.19  warships 
— /0.16  containerships 

H.  Bend.  Mom. 
(*a) 

0.85/0.15 

— 

Note:  In  ISSC  study  Xi,  Xi  are  normally  distributed  with  correlation  coefficient  0.7. 
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Table  4.7  Bounds  for  Bias  of  Response  Amplitude 
Operator  for  Hydrodynamic  Pressure  on 
SL-7  Containership  (Source:  Chen  et  al.  83) 


Froude  number 

Lower  Bound 

Upper  Bound 

0.15 

0.44 

1.35 

0.23 

0.41 

1.65 

0.32 

0.35 

1.60 

Using  Chen’s  results,  we  found  upper  and  lower  bounds  for  the  bias  in  the  response 
amplitude  operator  for  hydrodynamic  pressure.  (The  response  amplitude  operator  is  the 
square  root  of  the  ratio  of  the  spectral  ordinates  of  hydrodynamic  pressure  and  wave 
elevation  at  the  same  frequency.)  The  results  are  shown  in  Table  4.7  and  they  are  for  the 
SL-7  containeship. 

•  Clearly,  the  error  in  predicting  hydrodynamic  pressures  is  significantly  larger  than  that 
in  predicting  global  loads  (bending  moments  and  shear  forces).  This  is  true  because 
global  force  are  obtained  by  integrating  pressures  over  the  hull.  A  large  portion  of  the 
error  is  averaged  out  when  integrating.  Thus,  the  error  in  global  forces  is  smaller  than 
that  in  pressures. 

4.2  Offshore  platforms 

In  this  section,  we  study  uncertainties  in  short  term  and  long  term  extreme  loads  on 
offshore  platforms.  The  study  focuses  on  global  loads,  i.e.  base  shear  forces  and  overturning 
moments. 

Uncertainties  in  analysis  of  fixed  offshore  platforms  were  studied  in  the  context  of  the 
L.R.F.D.  A.P.I.  code  (Moses,  86).  The  maximum  annual  wave  height  was  assumed  lognor- 
mally  distributed.  Its  coefficient  of  variation  ranges  between  10  -  15%  for  the  North  Sea, 
15%  -  25%  for  the  Gulf  of  Mexico  and  it  is  somewhat  higher  for  offshore  Alaska  and  Cal¬ 
ifornia.  This  information  was  extracted  from  measurements  reported  by  various  authors 
and  it  is  summarized  on  p.  2.23  of  Moses  (86).  It  was  shown  that  the  effect  of  the  length 
of  the  exposure  time  on  the  lifetime  maximum  load  or  load  effect  is  to  reduce  its  COV 
and  to  increase  bias.  However,  no  information  was  provided  on  natural  and  subjective 
uncertainties.  Moreover,  these  two  types  of  uncertainty  were  not  distinguished  when  the 
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lifetime  distribution  of  the  maximum  wave  height  was  derived  from  the  distribution  of  the 
annual  maximum  wave  height. 

The  lifetime  maximum  platform  forces  were  assumed  to  be  related  to  the  maximum 
wave  height  according  to  the  following  relation, 

Fn  =  AH%  (4.10) 

where  A  is  called  analysis  coefficient,  and  the  exponent  a  is  1  and  2  for  inertia  and  drag 
dominated  platforms  respectively.  For  N  =  20  years,  the  bias  and  the  COV  of  the  analysis 
coefficient  were  assumed  to  be  0.93  and  0.25  respectively.  These  results  were  based  on 
measurements,  which  were  obtained  from  the  Ocean  Test  Structure  (O.T.S.)  (Anderson  et 
al  82).  The  latter  is  a  drag  dominated  platform. 

Olufsen  and  Bea  (90),  investigated  the  uncertainties  in  extreme  shear  force  and  over¬ 
turning  moment  for  two  platforms  located  in  the  Gulf  of  Mexico  and  in  the  North  Sea 
respectively.  They  considered  the  effects  of  both  random  (type  I),  and  modeling  uncer¬ 
tainties  (type  II).  An  empirical  model,  which  as  obtained  by  regression,  was  used  to  derive 
global  forces  from  the  wave  elevation.  The  following  uncertainties  were  taken  into  account: 

•  errors  in  the  procedure  for  deriving  the  force  from  the  extreme  wave  height, 

•  error  in  predicting  the  extreme  wave  height,  and 

•  uncertainties  due  to  the  effect  of  marine  fouling. 

It  was  stressed  by  Olufsen  and  Bea  (90)  that  the  coefficient  of  variation  in  the  extreme 
global  forces  is  severely  underestimated  if  modeling  uncertainties  are  not  treated  properly. 
More  specifically,  modeling  uncertainties,  which  are  involved  in  calculating  loads  and  their 
effects,  are  almost  perfectly  correlated  from  one  load  application  to  another.  Therefore,  in 
contrast  with  the  random  uncertainties,  modeling  uncertainties  do  not  decrease  with  the 
length  of  the  return  period  increasing.  Hence,  if  we  do  not  recognize  the  difference  between 
the  ways  that  the  above  two  types  of  uncertainties  propagate,  we  will  underestimate  the 
coefficients  of  variation  of  the  extreme  lifetime  loads.  It  is  striking  that  the  new  L.R.F.D. 
A.P.I.  design  code  is  based  on  a  value  of  0.37  for  the  coefficient  of  variation  of  the  20  year 
extreme  response  of  a  platform,  which  is  less  than  one  half  of  the  corresponding  value 
reported  by  Olufsen  and  Bea  (90)  (0.73  -  0.98  for  the  Gulf  of  Mexico,  and  0.65  for  the 
North  Sea). 

Bea  (89)  also  studied  uncertainties  for  a  platform  located  in  the  Northwest  Shelf  of 
Western  Australia.  The  global  force  to  the  platform  F,  was  calculated  by  the  following 
formula, 

F  =  kdkuHa  (4.11) 

where  kd  is  the  coefficient  in  the  relation  between  the  kinematics  of  the  water  particles  and 
F,  and  ku  denotes  the  coefficient  in  the  relation  between  the  former  and  the  wave  height, 


31 


H.  Thus,  the  product  kukd  corresponds  to  the  analysis  coefficient  in  eq.  (4.10)  which  was 
used  by  Moses  (86).  Exponent  a  is  1  and  2  for  inertia  and  drag  dominated  platforms, 
respectively.  Only  uncertainties  in  annual  maximum  values  were  reported.  These  values 
are  presented  in  Table  4.8.  The  bias  and  the  coefficient  of  variation, which  were  obtained 
by  combining  the  uncertainties  in  the  quantities  involved  in  calculating  global  forces  (eq. 
(4.11)),  were  found  to  be  in  good  agreement  with  the  corresponding  values  estimated  by 
direct  comparison  between  measurements  and  theoretical  predictions. 

The  principal  component  of  the  uncertainty  in  kd,  as  reported  by  Bea  (89),  is  uncertainty 
in  the  drag  coefficient  Cd  in  Morison’s  formula.  Based  on  OTS  data,  Bea  reported  that 
the  coefficient  of  variation  of  Kd  and  Cd,  which  corresponds  to  random  uncertainties  is 
0.10  and  that  due  to  modeling  uncertainties  are  0.23.  These  values  incorporate  the  effect 
of  marine  fouling. 

We  derived  the  uncertainties  in  the  lifetime  maximum  global  forces  from  those  of  the 
annual  maximum  loads  by  using  two  approaches,  in  order  to  demonstrate  how  important 
it  is  to  treat  modeling  uncertainties  properly.  In  the  first  approach,  (approach  1,  Table 
4.9)  we  assumed  that  modeling  uncertainties  are  perfectly  correlated  from  one  year  to 
another  and  we  used  equation  (2.2)  for  the  total  uncertainty.  In  the  second  approach, 
(approach  2,  Table  4.9)  we  treated  modeling  uncertainties  as  independent  from  one  year 
to  another.  The  details  on  the  calculation  of  the  coefficient  of  variation  are  described  in 
Appendix  A.  Return  periods  from  10  to  100  years  were  considered.  We  assumed  that 
random  variable  B/,  representing  random  uncertainties,  follows  the  lognormal  probability 
distribution.  Thus,  the  lifetime  maximum  value  of  Bj  follows  the  asymptotic,  type  /, 
extreme  probability  distribution.  (Appendix  A).  The  results  are  shown  in  Table  4.9.  It  is 
observed  that  the  second  approach  yields  significantly  lower  estimates  for  the  coefficient  of 
variation  than  the  first  approach.  This  is  because,  in  this  approach,  modeling  uncertainties 
are  treated  as  random  and,  therefore,  the  coefficient  of  variation  of  the  latter  decreases  with 
N  increasing.  Furthermore,  the  component  of  the  total  uncertainties  due  to  modeling  error 
does  not  change  with  N  in  the  first  approach.  It  is  observed  that  the  coefficient  of  variation 
of  the  maximum  force  over  a  20  year  period  is  0.66,  which  is  significantly  higher  than  the 
value  which  was  used  by  the  A.P.I.  rules  (0.37).  According  to  the  foregoing  discussion, 
this  discrepancy  is  due  to  the  difference  between  the  ways  that  modeling  uncertainties 
are  treated  in  this  study  and  in  Moses  (86).  However,  it  should  be  noted  that  the  actual 
numerical  values  for  extreme  load  might  be  different  than  the  values  reported  in  Table 
4.9,  because  it  is  difficult  to  distinguish  between  random  and  modeling  uncertainties,  and 
to  estimate  the  coefficients  of  variation  for  random  variables  B/  and  B/j.  However,  the 
trends  observed  in  Table  4.9  should  be  correct,  and  approach  1  is  more  appropriate  than 
approach  2  for  calculating  uncertainties  in  extreme  loads. 
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Table  4.8:  Uncertainties  in  Annual  Maximum  Loads 
for  Drag  and  Inertia  Dominated  Platforms  (Bea  89) 


Platform 

Quantity 

Random 

Modeling 

Type 

(type  I) 

(type  II) 

EBi 

COVB , 

EBn 

COVbu 

Wave  Height 

(H) 

1.0 

0.30 

1.1 

0.13 

Drag 

Domi¬ 

Kinematics 

(ku) 

1.0 

0.10 

0.41 

0.47 

nated 

Force  Coef. 

(*-) 

1.0 

0.10 

1.67 

0.23 

Global  Force 

1.0 

0.62 

0.83 

0.58 

Wave  Height 

(H) 

1.0 

0.30 

1.1 

0.13 

Inertia 

Domi¬ 

nated 

Kinematics 

(fcu) 

1.0 

0.10 

0.93 

0.20 

Force  Coef. 

(kd) 

1.0 

0.10 

0.65 

0.3 

Global  Force 

1.0 

0.33 

0.66 

0.38 

Table  4.9:  Total  Coefficient  of  Variation  of  Global  Force 
as  a  Function  of  Return  Period 


Years  Approach  1  Approach  2 

COVB,  COVb  COVb 


Init.  Dist.  (Corr.  to  1  yr.) 

0.62 

0.85 

0.85 

10 

0.35 

0.68 

0.43 

20 

0.31 

0.66 

0.38 

50 

0.27 

0.64 

0.34 

100 

0.25 

0.63 

0.31 

Wirsching  (81)  also  studied  uncertainties  in  loads  applied  to  offshore  platforms.  He 
represented  the  uncertainties  in  loads  by  the  product  of  two  coefficients  denoted  by  Bs 
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and  Bf •  Bs  corresponds  to  environmental  uncertainties  and  Bf  accounts  for  the  error  in 
load  calculation.  The  statistics  of  Bs  were  presented  in  chapter  3.  Bf  was  assumed  to  be 
lognormally  distributed  with  an  average  value  ranging  from  0.6  to  1.1  and  a  COV  between 
0.1  and  0.3. 

Guedes  Soares  and  Moan  (82)  considered  the  uncertainties  in  the  extreme  forces  applied 
to  a  vertical  pile  in  the  North  Sea.  The  extreme  forces  correspond  to  a  return  period  of 
100  years.  Table  4.10  presents  the  random  variables  which  were  considered  in  this  study 
and  their  means  and  COV’s. 

The  COV  of  the  extreme  load  was  found  to  be  in  the  range  between  0.34  and  0.45.  The 
uncertainty  in  the  extreme  wave  height  was  found  to  be  the  most  important,  because  its 
effect  on  the  global  load  was  considerably  larger  than  the  effects  of  all  the  other  uncertain¬ 
ties.  This  conclusion  agrees  with  the  conclusions  from  Wirsching  and  Bea.  Therefore,  the 
uncertainty  in  environmental  description  (the  extreme  wave  height)  is  the  most  important 
for  offshore  platforms. 

The  modeling  uncertainties  in  estimating  wind,  current  and  tide  loads  on  tension  leg 
platforms  were  quantified  in  the  context  of  a  study  which  was  sponsored  by  ABS  and 
Conoco  (Wirsching  82).  The  bias  associated  with  the  effect  of  wind,  current  and  tide  loads 
was  assumed  to  be  the  product  of  five  random  variables  which  account  for  uncertainties  in 
predicting  the  most  critical  combination  of  wind  current  and  tide  conditions,  calculation  of 
induced  loads,  calculation  of  dynamic  response,  calculation  of  nominal  member  loads,  and 
structural  analysis.  These  r.v.’s  were  assumed  to  be  lognormally  distributed.  The  median 
biases  and  COV’s  of  those  random  variables  are  presented  in  Table  4.11. 

The  values  in  Table  4.11  are  based  on  subjective  estimates  of  experts.  It  is  mentioned 
that  the  estimates  for  the  COV’s  might  be  low  (Wirshing  82).  The  correlation  between 
the  bias  in  wind,  current  and  tide  and  in  wave  loads  was  not  estimated  in  the  above  study 
although  it  is  recognized  that  it  might  be  significant  in  developing  probabilistic  design 
rules. 

The  results  from  the  studies  considered  in  this  section  are  summarized  in  Table  4.12. 
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Table  4.10  Uncertainties  Involved  in  Predicting  Extreme 
Loads  on  a  Vertical  Pile  in  the  North  Sea 
(Source:  Guedes  Soares  and  Moan  (82)) 


Quantity 

Mean  value 

COV 

Extreme  wave  height  (H) 

30m 

0.16 

Wave  period  (T) 

5.4  +  0.373  H  (sec) 

0.14 

Water  depth  ( D ) 

80m 

2/D 

Current  velocity  (C) 

1.25  m/sec 

0.35 

Pile  diameter  ( D ') 

4.0m 

0.0 

Fouling  thickness  ( K ) 

0.175m 

0.45 

Surface  roughness  (R) 

0.02 

0.4 

Drag  coefficient  (Cd) 

Sarpkaya’s  data 

0.1 

Inertia  coefficient  (Cm) 

Sarpkaya’s  data 

0.1 

Wave  Kinematics 

Stokes  theory 

0.25 

Note:  The  following  correlation  coefficients  were  assumed  for  the  above  random  variables: 
P(H,T)=  0.5, 
p(//,C)=0.4, 
p(K,R)  =  0.7, 
p(R,Cd)=  0.5, 

P  (R,Cm)  =  -0.5, 

P  ( Cd ,  Cm)  =  -0.9. 
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Table  4.11  Uncertainties  in  Wind,  Current  and 
Tide  Loads  (Wirshing  82) 


Type  of  Uncertainty 

Bias  (Median) 

cov 

Prediction  of  most  critical  combination 
of  wind,  current  and  tide  loads 

0.95 

0.12 

Calculation  of  induced  load 

0.92 

0.10 

Calculation  of  dynamic  response 

0.95 

0.14 

Calculation  of  nominal  number  loads 

0.90 

0.15 

Structural  analysis 

1.0 

0.05 

Table  4.12  Uncertainties  in  Extreme  Global  Loads 
on  Offshore  Platforms 

Source  Bias/COV  Return  Period 

(years) 

20 

20 
100 

100 

100 
100 

Notes: 

1  Bias  is  lognormally  distributed 

2  Gulf  of  Mexico 

3  North  Sea 


Moses  (86) 
Bea  (89) 


Wirsching  (81) 
Guedes  Soares  (82) 
Olufsen  and  Bea  (90) 


0.7/0.37 

-/0.661 
— /0.631 

0.4  -  1.3/0.41  -  0.671 

— /0.34  -  0.45 

— /0.73  -  0.931,2 
-/0.651’3 
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5.0  Load  Combinations 

In  many  cases,  two  or  more  loads  or  load  effects  act  on  the  same  component  of  a  marine 
structure.  It  is  very  important  to  estimate  the  extreme  values  of  the  combined  load  in 
order  to  evaluate  the  safety  of  the  structure  and  make  design  decisions.  The  problem  of 
combining  loads  and  of  finding  the  extreme  values  of  the  combined  load  might  be  very 
difficult  because,  for  most  of  the  cases,  the  loads,  which  must  be  combined,  vary  randomly 
in  time  and  their  maximum  values  do  not  occur  at  the  same  time  instant.  Moreover,  it  is 
difficult  to  determine  the  location  in  the  structure  where  the  combined  load  effect  takes 
its  maximum  value.  In  this  chapter,  we  investigate  the  errors  associated  with  combining 
loads  in  ships. 

Examples  of  loads  effects,  which  need  to  be  combined  are  the  following: 

a)  Stillwater  and  vertical  wave  bending  moments, 

b)  vertical  and  horizontal  wave  bending  moments, 

c)  wave  bending  and  torsional  moments, 

d)  wave  and  slamming  bending  moments. 

Guedes  Soares  (84)  has  shown  that  the  Ferry  Borges  and  Castanheta  (68),  and  Larrabee 
and  Cornell  (81)  methods  are  almost  exact  for  problem  a).  Problems  b)  and  c)  can  be 
easily  solved,  provided  that  the  cross  spectral  densities  of  the  pertinent  random  processes 
are  known.  Although  this  can  be  easily  done  by  using  any  linear  seakeeping  computer 
code,  most  of  the  studies  on  longitudinal  ship  strength  do  not  address  the  issue  of  load 
combination. 

Problem  d)  is  more  difficult  because  the  load  processes,  which  must  be  combined,  are 
correlated,  and  the  slamming  bending  moment  is  nonstationary.  Moreover,  it  is  very 
difficult  to  calculate  the  relative  positions  of  the  wave  and  slamming  bending  moments 
peaks,  in  the  time  scale. 

In  order  to  assess  the  uncertainties  in  combining  slamming  and  wave  bending  moments, 
we  followed  a  Monte  Carlo  simulation  approach.  More  specifically,  we  developed  a  com¬ 
puter  code,  which  calculates  slamming  and  wave  bending  moments  and  combines  them  for 
a  given  time  history  of  the  wave  elevation.  The  code  is  an  implementation  of  the  procedure 
suggested  by  Kaplan  (72,87)  for  generating  random  time  histories  of  the  wave  elevation, 
and  the  method  described  by  Kaplan  (72,86)  for  calculating  the  resulting  slamming  bend¬ 
ing  moments.  We  generated  a  large  number  of  wave  elevation  time  histories,  and  for  each 
of  them,  we  calculated,  a)  the  time  histories  of  the  resulting  slamming,  wave,  and  combined 
bending  moments,  and  b)  the  maxima  of  the  combined  bending  moments.  For  the  purpose 
of  this  study,  we  assumed  that  these  are  the  actual  values  of  the  maxima  of  the  combined 
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bending  moment.  We  also  estimated  the  maxima  of  the  combined  bending  moment,  from 
those  of  the  slamming  and  wave  bending  moments,  by  using  standard  methods.  These 
methods  are,  Turkstra’s  rule,  the  peak  coincidence  approximation,  and  the  square  root  of 
sum  of  squares  (S.R.S.S.)  rule,  and  they  are  described  in  Appendix  B.  Then,  we  compared 
the  actual  maxima  against  those  from  the  aforementioned  methods.  The  results  of  this 
study  are  summarized  in  Table  5.1.  The  details  of  our  approach  can  be  found  in  Appendix 
B. 


The  following  conclusions  can  be  extracted  from  Table  5.1,  and  from  the  results  pre¬ 
sented  in  Appendix  B: 

a)  Turkstra’s  rule  is  unconservative.  For  some  cases,  it  may  underestimate  the  combined 
bending  moment  by  40%. 

b)  The  peak  coincidence  approximation  is  overconservative,  and  as  such,  it  might  lead  to 
very  inefficient  designs. 

c)  The  bias  of  the  S.R.S.S.  approximation  was  found  to  be  almost  1.00.  However,  this 
approximation  has  slightly  larger  scatter  than  the  other  two  methods.  Moreover,  it  is 
conservative  for  some  cases,  while  it  is  unconservative  for  some  others. 

d)  The  uncertainty  associated  with  the  combination  of  the  above  moments  is  significant. 

e)  Presently,  there  is  no  method,  which  can  be  used  as  a  reliable  engineering  tool  to  predict 
the  maxima  of  the  combined  slamming  and  wave  bending  moments. 

Table  5.1  Bias  and  COV  of  Standard  Methods  for 
Combining  Slamming  and  Wave  Bending  Moments 


Method 


Bias 


COV 


Turkstra’s  rule 

1.17 

0.11 

Peak  coincidence 

0.72 

0.11 

approximation 

SRSS 

1.01 

0.12 

38 


6.0  Structural  analysis 

In  this  chapter  we  investigate  the  accuracy  and  the  fidelity  of  procedures  which  are 
used  to  predict  the  response  of  marine  structures  to  the  applied  loads.  We  will  focus 
on  uncertainties  associated  with  finite  element  analysis  (F.E.A.)  procedures,  but  we  will 
also  examine  the  accuracy  of  traditional  beam  theory  in  estimating  main  girder  stresses 
in  the  ship  hull.  Information  is  also  provided  on  uncertainties  associated  with  F.E.A.  of 
automotive  and  aerospace  structures. 

6.1  Ships 

6.1.1  Uncertainties  in  F.E.A. 

The  accuracy  of  F.E.A.  procedures  can  be  assessed  by  comparing  their  results  against 
measurements  from  test  models  or  actual  ships.  Elbatouti  et  al  (76),  Webster  and  Payer 
(77),  and  Westin  (81)  presented  experimental  and  analytical  results  for  the  SL-7  containe- 
ship.  Analytical  results  were  obtained  by  F.E.A.,  or  simplified  approaches  such  as  the 
Finite  Beam  Technique  (Westin  81)  and  the  beam  theory.  Experimental  measurements 
were  obtained  from  tests  which  were  carried  out  on  a  steel  structural  model.  Jan  et  al 
(79)  also  presented  experimental  and  analytical  results  on  wave  bending  stresses  in  the 
hull  of  an  SL-7  containership.  Here,  some  of  the  above  results  are  analyzed  to  assess  the 
uncertainties  in  structural  analysis. 

A  sketch  of  the  ship  hull  including  the  locations,  at  which  measurements  were  taken  in 
Elbatouti’s  study,  is  presented  in  Figure  6.1.  We  considered  two  loading  cases,  which  are 
shown  in  Figure  6.2.  In  the  first  case,  the  ship  was  subjected  to  a  vertical  bending  moment 
(Loading  case  1).  In  the  second  case,  the  ship  was  subjected  to  pure  torsion.  (Loading 
case  2).  A  finite  element  model  F.E.M.  of  the  SL-7  ship  hull  was  generated  and  analyzed 
by  Elbatouti  by  using  the  ABS/DAISY  computer  code.  This  model  consisted  of  roughly 
3,000  elements,  which  were  mostly  triangular  plates.  Transverse  bulkheads  were  modelled 
by  membrane  and  bar  elements.  The  model  had  1800  nodes  and  6700  degrees  of  freedom. 

Figure  6.3  compares  measurements  and  predictions  of  the  angle  of  twist  along  the  hull 
under  loading  case  2.  It  can  be  observed  that  the  F.E.M.  is  stiffer  than  the  steel  model. 
This  is  a  general  conclusion.  For  the  majority  of  the  cases  in  all  engineering  applications, 
F.E.A.  underestimates  flexibility.  It  can  also  be  observed  that  there  is  a  sudden  change  in 
the  slope  of  the  twisting  angle  at  point  A.  (This  point  corresponds  to  frame  290).  This 
change  is  due  to  the  large  change  in  the  torsional  stiffness  of  the  hull  at  frame  290.  The 
discrepancy  between  the  slopes  of  the  measured  and  experimental  twisting  angles  is  large 
at  this  frame.  This  is  probably  due  to  the  fact  that  hull  warping  is  not  constrained  in  the 
F.E.M.  model  as  strongly  as  in  the  steel  model.  Based  on  the  results  in  Figure  6.3,  the 
bias  in  the  twist  angle  was  found  equal  to  1.4  and  the  COV  equal  to  0.18. 
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Figures  6.4  and  6.5  depict  the  normal  warping  stresses  in  frames  142  and  290,  respec¬ 
tively  under  loading  case  2.  Analytically  predicted  stresses,  as  well  as  those  measured  from 
tests  on  a  steel  and  a  vinyl  model,  are  plotted  in  these  Figures.  The  bias  for  frame  142  was 
estimated  1.20  and  the  COV  0.42.  These  estimates  are  based  on  the  analytical  stresses 
and  those  obtained  from  tests  on  the  steel  model.  For  frame  290,  these  numbers  were  1.8 
and  0.67,  which  indicates  that  the  error  in  the  warping  stresses  was  larger  at  this  frame. 
This  is  attributed  to  the  sudden  change  in  the  geometry  and  torsional  stiffness  of  the  hull 
at  this  frame.  Moreover,  according  to  the  Figure  6.3  and  the  discussion  in  the  previous 
paragraph,  the  F.E.M.  does  not  represent  accurately  the  warping  constraint  at  frame  290. 

Clearly,  it  is  very  difficult  to  predict  the  normal  warping  stresses  in  ships  with  large 
openings  in  their  decks,  by  using  F.E.A.  Even  if  we  use  detailed  F.E.M’s  to  represent 
portions  of  the  structure,  at  which  the  geometry  suddenly  changes,  the  error  is  still  likely 
to  be  large.  The  reason  is  that  we  cannot  accurately  determine  the  forces  applied  to  the 
boundaries  of  these  critical  portions  of  the  hull  from  the  rest  of  the  ship.  Analysis  of  frame 
290  in  the  SL-7  containeship,  is  an  example  of  a  case,  for  which  F.E.A.  cannot  accurately 
predict  the  normal  forces  on  this  frame  due  to  the  effect  of  constraining  warping.  (Elbatouti 
76). 

The  shear  stress  distribution  at  frame  186  under  loading  case  2  is  plotted  in  Figure 
6.6.  We  observe  that  the  F.E.A.  results  are  in  good  agreement  with  measurements.  The 
bias  was  found  equal  to  1.01,  which  indicates  that  there  is  no  systematic  trend  for  the 
predictions  to  be  either  conservative  or  unconservative.  The  COV  was  found  0.20,  which 
implies  that  shear  stress  predictions  are  more  accurate  than  those  for  warping  stresses. 

Finally,  we  compared  predicted  and  measured  normal  longitudinal  stresses  under  pure 
bending  (loading  case  1).  The  results  are  plotted  in  Figure  6.7  together  with  those  from 
beam  theory.  The  bias  and  the  COV  for  F.E.A.  results  are  0.93  and  0.17,  respectively. 
This  implies  that  results  on  longitudinal  bending  stresses  are  more  accurate  than  those 
for  shear  or  warping  stresses.  Moreover,  beam  theory  was  found  to  be  very  accurate  in 
predicting  normal  stresses.  Its  bias  was  0.94  and  its  COV  was  0.10. 

It  should  be  emphasized  that  tb  zesults  of  the  previous  paragraphs  are  based  on  a 
small  number  of  samples.  Therefore,  they  should  be  used  with  caution.  Furthermore,  the 
SL-7  is  a  very  difficult  ship  to  analyze  due  to  the  extremely  complex  geometry  of  its  hull. 
The  COV’s  for  the  results  from  F.E.A. ’s  of  other  ships,  such  as  tankers,  are  expected  to 
be  considerably  lower  than  the  values,  which  were  estimated  in  this  example.  Moreover, 
the  data,  which  were  used  to  estimate  uncertainties,  were  contaminated  with  experimental 
error.  Unfortunately,  we  cannot  get  rid  of  this  type  of  error.  Thus,  the  values  of  the 
COV’s  are  likely  to  be  conservative  because,  they  incorporate  the  effects  of  the  error  in 
both  F.E.A.  and  experimental  measurements. 

Jan  et  al  (79)  compared  theoretically  predicted  and  measured  stresses  in  the  hull  of 
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Figure  6.2  Loading  Conditions 
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Figure  6.3  Comparison  of  Twisting  Angles  Along 
the  Ship  Hull 
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Figure  6.4  Longitudinal  Warping  Stresses  in 
Frame  142 
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Figure  6.6  Shear  Stresses  at  Frame  186 
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Figure  6.7  Normal  Bending  Stress  in  Frame  186 
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an  SL-7  containership.  Strip  theory  based  computer  codes  (ABS/SHIPMOTION  and 
ABS/DUNPRE)  were  used  to  predict  hydrodynamic  motions  and  loads.  The  ABS/DAISY 
code  was  used  for  F.E.A.  Stress  spectra  and  rms  of  wave  bending  stress  in  head  and  obligue 
seas  were  investigated.  The  predicted  rms  of  the  vertical  bending  stress  in  head  seas  was 
also  compared  against  measurements.  The  stress  was  calculated  from  the  spectrum  of  the 
wave  elevation  and  the  stress  response  amplitude  operator,  which  was  obtained  analytically. 
The  results,  as  well  as  the  estimated  bias  for  the  rms  of  the  bending  stress  amidships  are 
shown  in  Table  6.1. 


Table  6.1  Measured  and  Predicted  Longitudinal 
Wave  Bending  Stresses  Amidships 
(Source:  Jan  et  al,  79) 


Calculated 

(RMS) 

Measured 

(RMS) 

B 

(RMS) 

Comments 

8343 

6743 

0.81 

the  results  in 

7295 

6344 

0.87 

columns  1  and  2 

7183 

5368 

0.75 

refer  to  the 

3547 

2990 

0.84 

RMS  of  peak 

4192 

3300 

0.79 

to  trough 

4544 

4730 

1.04 

stresses  on 

2542 

2000 

0.79 

main  deck 

4252 

3410 

0.80 

amidships 

4766 

4310 

0.90 

3188 

3530 

1.11 

4595 

5800 

1.26 

Based  on  the  results  in  Table  6.1,  the  bias  was  found  equal  to  0.91  and  the  COV  0.16. 
These  values  incorporate  errors  in  both  the  response  amplitude  operator  between  the  wave 
elevation  and  wave  loads,  and  the  structural  analysis  of  the  hull.  Assuming  that  the  COV 
for  the  response  amptitude  operator  between  the  wave  elevation  and  the  bending  moment 
is  0.1  (Kaplan  84),  the  COV  for  the  stress  analysis  is  estimated  0.125.  The  COV  is  of  the 
same  order  of  magnitude  with  that  found  from  analysis  of  Elbatouti’s  results. 

The  results  on  uncertainties  in  structural  analysis  of  ship  hulls  are  summarized  in  Table 

6.2. 
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Table  6.2  Uncertainties  in  F.E.A.  of  an  SL-7  Containeship 


Loading  Case 

Quantity 

Bias 

COV 

Pure  longitudinal 
bending 

Normal  stress 

0.93 

0.12  -  0.17 

Distributed  tor¬ 
sional  moment 

Shear  stress 

1.01 

0.20 

Distributed  tor¬ 
sional  moment 

Normal  warping 
stresses 

1.20 

0.42 

Distributed  tor¬ 
sional  moment 

Twisting  angle 
along  the  hull 

1.4 

0.18 

6.1.2  Uncertainties  in  F.E.A.  in  other  engineering  applications 

It  is  interesting  to  compare  the  errors  in  Table  6.2  with  those  associated  with  F.E.A. 
of  other  engineering  structures.  Sutharsana  et  al  (90)  and  Newlin  et  al  (90)  reviewed  the 
literature  on  F.E.A.  of  aerospace  structures  and  discussed  the  subject  with  experts.  Their 
conclusion  was  that  F.E.  stress  analysis  would  calculate  stresses  to  within  20%  of  the  true 
value.  Based  on  this  conclusion,  they  assumed  that  the  analysis  bias  follows  the  uniform 
probability  distribution  from  0.8  to  1.2.  According  to  this  assumption,  the  average  bias  is 
1.0  and  the  COV  is  0.12. 

F.E.A. ’s  using  F.E.M.  with  roughly  5000  degrees  of  freedom  are  routinely  used  to  ana¬ 
lyze  automotive  structures.  The  author  had  extensive  discussions  with  experts  on  F.E.A. 
of  car  bodies  on  the  correlation  between  analytical  predictions  and  measurements.  The 
following  conclusions  are  based  on  these  discussions:  It  is  general  consensus  that  F.E.A. 
underestimates  the  flexibility  of  a  car  body.  This  conclusion  agrees  with  that  for  F.E.A. 
of  ship  hulls  (Table  6.2).  Moreover,  the  error  in  predicting  deflections  due  to  bending 
or  torsional  loads  ranges  between  10  to  20%.  Finally,  it  is  more  difficult  to  predict  the 
response  due  to  longitudinal  bending  than  that  due  to  torsional  loading. 

6.1.3  Errors  Due  to  Shear  Lag 

Elementary  theory  of  bending  assumes  that  the  bending  stresses  in  a  beam  are  pro¬ 
portional  to  the  distance  from  the  neutral  axis.  This  is  not  true  for  sections  having  wide 
flanges.  Indeed,  the  parts  of  the  flange  fax  away  from  the  web  can  carry  lower  stresses, 
which  makes  the  actual  maximum  normal  stress  larger  than  that  predicted  by  the  classical 
beam  theory. 


For  design  purposes,  this  is  dealt  with  by  reducing  the  width  of  the  flange  so  that  the 
correct  bending  stress  can  be  obtained  by  applying  the  elementary  theory  of  bending. 


In  the  literature,  the  terms,  “effective  breadth”  and  “effective  width”  are  used  inter¬ 
changeably  for  two  different  phenomena: 

(i)  Effective  width  of  plates  subjected  to  compressive  loads  to  compensate  for  buckling 
effects. 


(ii)  Effective  width  of  plates  subjected  to  lateral  loads  to  compensate  for  shear  lag  effects. 

Schade  (51,53)  attempted  to  resolve  this  confusion  by  using  the  terms  “effective  width” 
for  case  (i)  and  “effective  breadth”  for  case  (ii). 

Here,  we  review  various  methods  for  estimating  effective  breadth.  The  methods  are 
applied  to  some  simple  plate-stiffener  assemblies  subjected  to  bending,  and  their  results 
are  compared  to  test  measurements  in  order  to  quantify  the  bias  of  each  method.  The 
following  methods  are  studied: 

a.  Timoshenko  fe  Goodier 


A  stress  function  is  defined  such  that, 


d2<j> 

d2(f> 

ay~lh? 

d2<t> 

T~  dxd. 

The  stress-strain  relations  and  the  compatibility  condition  together  result  in  the  bihar¬ 
monic  equation 

gv  ,  g  0V  I  ay  n 

dx 4  dx2dy2  dy 4 

An  expression  for  4>(n ,  y)  satisifying  the  boundary  conditions  is  obtained  from  this  equation 
from  which  the  stresses  are  calculated.  The  effective  breadth  is  then  calculated  by 


A 


where  as  is  the  stress  in  the  part  of  the  flange  close  to  the  web.  In  case  of  a  general 
transverse  loading,  the  bending  moment  is  represented  by  a  sine  series  of  the  type: 

M(x)  =  EM„sin  nnx/L 
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where 


Mn  =  -  f  M (x)sm^^-dx 
wo  L 

This  yields  the  following  formula  for  effective  breadth: 

M{x) 


7!  Mn  sinfrura/L) 
n=l,3,5  4+(k/ P)nx 


—  4]  where, 


k  =  i lUH&zul 

4  * 

Aa  is  the  stiffer  area. 
t  is  the  plate  thickness, 

L  is  the  beam  length,  and 
H  is  Poisson’s  ratio, 
b.  Schade 

According  to  Schade,  the  effective  breadth  depends  upon  the  following  factors: 

(a)  Span  of  the  stiffners.  It  increases  with  the  span  increasing. 

(b)  Nature  of  the  load,  but  not  its  magnitude.  It  is  lowest  where  high  shear  exists. 

(c)  The  boundary  conditions,  particularly  at  the  plate  sides.  It  does  not  depend  on 
thickness. 

For  a  multiple  stiffener  configuration, 


A 

b 


1.1 


where  CL  ~  distance  between  points  of  zero  bending  moment. 


c.  Bureau  of  Ships  (BuShips) 

This  study  assumes  that  the  effective  breadth  is  influenced  not  only  by  the  thickness 
of  the  plating,  the  spacing,  and  the  span  of  the  stiffners  but  also  by  the  shape  factor  of 
the  stiffner  as  well  as  the  general  stress  level.  This  approach  has  been  presented  in  Schade 
(51)  (pp.  421-422). 


44 


The  equation  which  relates  these  terms  is, 

a  IT  1, _ sinhf? 

9  \  —  _  y  ^  _  a  ] 

0yp  Cc  et  sinh^(l  +  £  +  ^^r{(3  -  p)(l  +  fx)cosh^  -  (1  + /i)2$cosh$} 
where, 

t  -  thickness  of  plate, 

p  -  Poisson’s  ratio, 

2b  -  breadth  of  plate  or  flange, 

2A  -  effective  breadth, 

a  -  ultimate  compressive  stress  of  panel, 

cryp  -  yield  stress, 

It  -  Total  Moment  of  Inertia, 

Cc  -  distance  of  neutral  axis  from  top  fiber  of  plate, 

e  -  distance  from  N.A.  of  stiffner  alone  to  c.g.  of  plate, 

p  -  radius  of  gyration  of  stiffner, 

Is  -  M.I.  of  stiffner,  and 

2 V  —  L  —  Span  of  the  Beam. 

d.  Modified  Schade's  Formula  (Faulkner  75) 

Schade’s  formula  was  modified  as  follows: 

A  1.1c 
*"1  +  2®  2 

where  c  is  constant  equal  to  0.75. 

In  the  following  we  estimate  the  bias  and  scatter  introduced  by  the  aforementioned 
formulas. 

Consider  a  simply  supported  beam  under  a  uniformly  distributed  load  along  its  length. 
The  beam  consists  of  a  plate  and  a  stiffener,  which  is  attached  to  the  beam.  The  dimensions 
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are  provided  in  Tables  6.3  and  6.5.  A  comparison  of  the  previous  methods  and  some  test 
results  has  been  made  in  Tables  6.4,  6.6  and  6.7. 

Table  6.3  Loading  and  Geometry  Description: 

Loading  -  Uniformly  Distributed. 

Span  -  16’ 

Spacing  (Plate  breadth)  -  24” 

Stiffener  Area  -  6.22  sq.  in. 

The  following  table  presents  the  effective  breadth  for  various  plate  thickness. 

Table  6.4  Effective  Breadth  For  Various  Plate  Thicknesses 

Effective  Breadth 


Plate  Thickness 
(in) 

Tests 

Schade 

Modified 

Schade 

Timoshenko 

BuShips 

0.32 

0.60 

1.06 

0.80 

0.77 

0.62 

0.44 

0.88 

1.06 

0.80 

0.84 

0.89 

0.63 

0.89 

1.06 

0.80 

0.70 

0.87 

Test  results  are  taken  from  Schade  (51) 

Table  6.5  Loading  and  Geometry  Description: 

Loading  -  Uniformly  Distributed 
Span  -  16’ 

Spacing  -  24” 

Plate  Thickness  -  0.44” 

The  following  table  shows  the  effective  breadth  for  various  stiffeners  with  different  cross- 
sectional  areas. 


46 


Table  6.6:  Effective  Breadth  For  Various  Stiffeners 
with  Different  Cross-sectional  Areas 


Modified 


of  Stiffness 
in 

Tests 

Schade 

Schade 

Timoshenko 

BuShips 

6.22 

0.88 

1.06 

0.80 

0.84 

0.89 

7.72 

0.82 

1.06 

0.80 

0.99 

0.82 

8.61 

0.82 

1.06 

0.80 

1.06 

0.87 

10.72 

0.79 

1.06 

0.80 

1.16 

0.81 

The  aforementioned  results  were  analyzed  in  order  to  estimate  the  bias  and  scatter  of 
the  effective  breadth  and  the  maximum  stress  as  obtained  from  each  formula.  The  results 
are  presented  in  the  following  table. 

Table  6.7:  Bias  and  Scatter  of  Methods  to  Estimate  Effective  Breadth 


Method 

Stress 

Bias 

Scatter 

Schade 

1.35 

0.16 

Timoshenko 

1.17 

0.21 

BuShips 

1.02 

0.03 

Beam  Theory 

1.27 

0.16 

Modified  Schade 

1.02 

0.16 

Clearly,  the  effect  of  shear  lag  may  be  significant  for  some  cases.  This  is  particularly 
true  for  longitudinal  bending  of  short  beams  (length/breadth  <  5),  for  which  case  shear 
effects  are  important.  Transverse  bending  of  a  SWATH  ship  is  an  example  for  which  shear 
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effects  must  be  taken  into  account.  The  results  in  Table  6.7  demonstrate  that  classical 
beam  theory  underestimates  bending  stresses.  Empirical  equations,  such  as  the  modified 
Schade’s  formula,  can  be  used  to  account  for  shear  lag  in  stress  analysis.  Moreover,  one 
of  the  aforementioned  empirical  equations  (BuShips)  was  found  to  be  very  effective  in 
calculating  effective  breadth  and  the  maximum  stress. 

6.1.4  Effect  of  joint  flexibility  and  rigid  beam  lengths  in  transverse  strength 
analysis  of  ships. 

In  strength  analysis  of  traverse  hull  frames,  it  is  assumed  that  the  joints,  which  connect 
frame  elements,  are  rigid.  Moreover,  it  is  assumed  that  transverse  frames  are  flexible  all 
over  their  length.  However,  in  reality,  joints  are  flexible,  and  portions  of  transverse  frames 
in  the  vicinity  of  joints  are  practically  rigid  because  they  are  reinforced  by  brackets  and 
by  the  attached  longitudinal  beams. 

Hughes  (83)  studied  the  effect  of  joint  flexibility  and  rigid  beam  lengths.  His  study 
was  based  on  Milchert’s  (72)  results.  It  was  demonstrated  that  an  idealization  of  the 
transverse  frames  of  a  ship  as  continuous  frames  with  no  allowance  for  joint  flexibility 
and  rigid  frame-length  may  lead  to  significant  errors  in  transverse  strength  analysis.  The 
results  reported  by  Milchert  (72)  and  Hughes  (83)  showed  that  the  bending  moment  which 
is  calculated  by  the  aforementioned  idealization,  might  be  50%  larger  than  the  actual  one. 
Moreover,  the  displacements  might  be  overestimated  by  60%.  However,  the  significance  of 
these  results  is  reduced  by  the  fact  that  they  correspond  to  a  frame  whose  members  are 
very  short  (the  size  of  the  joint  was  roughly  30%  of  that  of  the  frame  element).  For  this 
case,  it  is  not  necessary  to  be  very  accurate  when  analyzing  transverse  strength,  because 
the  displacements  and  stresses  are  small.  Furthermore,  the  errors  due  to  joint  flexibility 
and  rigid  frame  lengths  should  be  considerably  smaller  than  for  frames  with  long  elements. 

The  effect  of  joint  flexibility  and  rigid  frame  lengths  can  be  accounted  for  if  we  employ  a 
very  detailed  finite  element  model,  which  models  transverse  frames  and  brackets  by  plate 
elements.  The  computer  program  MAESTRO  (Hughes  83),  which  is  for  ship  structural 
analysis,  accounts  for  the  joint  flexibility  and  the  effect  of  rigid  frame  lengths.  Most  of  the 
general  purpose  finite  element  codes  do  not  account  for  the  above  effects. 

6.2  Offshore  platforms. 

Finite  element  programs,  such  as  NASTRAN,  are  used  to  derive  member  forces  in 
offshore  platforms  from  the  applied  loads.  Wirsching  (81)  reported  that  the  bias,  which 
is  associated  with  the  calculation  of  member  forces,  ranges  between  0.8  and  1.10,  while 
the  COV  is  between  0.2  to  0.4.  These  numbers  were  based  on  discussions  with  experts  in 
structural  analysis  of  platforms. 

The  error  in  structural  analysis  was  attributed  to  the  effect  of  joint  flexibility,  which 
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is  not  taken  into  account  in  F.E.A.,  and  uncertainty  in  damping.  The  COV  of  member 
force  estimates  due  to  the  latter  effect  was  estimated  as  8%.  Due  to  the  joint  flexibility, 
the  actual  natural  period  of  a  platform  is  lower  than  that  estimated  by  F.E.A.  Thus, 
joint  flexibility  introduces  uncertainties  in  the  natural  period.  Moreover,  by  treating  joints 
as  rigid,  we  overestimate  bending  moments  and  the  resulting  stresses  in  the  structural 
members,  which  are  attached  to  the  joints.  This  effect  was  also  observed  in  the  analysis 
of  transverse  ship  frames  (Hughes,  83). 

Bouwkamp  et  al.  (80),  found  that  the  effect  of  joint  flexibility  on  the  structural  response 
of  an  offshore  platform  might  be  significant  for  high  offshore  towers  (height  >  500  ft). 
Moreover,  the  effect  of  joint  flexibility  on  member  forces  is  strong  for  tubular  members 
whose  lengths  are  short  compared  to  their  diameter.  It  is  interesting  that,  although  joint 
flexibility  increases  the  deflections  of  a  platform,  classical  structural  analysis  procedures, 
which  model  joints  as  rigid,  might  overestimate  the  displacements  at  some  locations.  The 
maximum  error  in  deflection  due  to  joint  flexibility  was  found  equal  to  30%  for  the  example 
considered  in  Bouwkamp  (80).  The  corresponding  maximum  error  in  the  bending  moments 
in  structural  members  was  found  to  range  up  to  500%.  The  error  in  the  first  natural  period 
was  small,  less  than  2%.  However,  the  error  in  higher  natural  periods  was  higher  (for 
example  it  is  11%  for  the  second  natural  period).  These  modest  errors  in  natural  period 
are  important  in  structural  analysis  because  they  lead  to  large  errors  in  the  dynamic 
response. 

Moses  (85),  investigated  uncertainties  in  dynamic  analysis  of  fixed  offshore  platforms 
with  low  natural  frequencies,  for  which  dynamic  effects  are  important.  He  used  the  inertia 
load  factor  (I.L.F.)  in  order  to  account  for  dynamic  effects.  The  I.L.F.  is  defined  by  the 
following  equation; 

Dynamic  response  =  (1.0  +  ILF)x  static  response. 

The  overall  uncertainty  of  ILF  was  estimated  by  using  a  FOSM  approximation  and  Taylor 
expansion  about  the  mean  values  of  the  random  variables,  which  are  used  to  quantify 
uncertainties  in  all  the  steps  of  dynamic  analysis.  The  COV  of  the  ILF  for  member  loads 
was  found  equal  to  0.61.  Table  6.8  summarizes  the  uncertainties  which  were  considered  by 
Moses. 
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Table  6.8:  Uncertainties  in  Dynamic  Analysis  of 
Fixed  Offshore  Platforms 


Quantity 

Uncertainty 

(COV) 

Effect  on  ILF 

Relative  contribution 
(%) 

Natural  period 

0.1 

0.48 

62 

Damping 

0.3 

0.10 

2 

Spectral  shape  variability 

- 

0.20 

11 

Transfer  function 

- 

0.20 

2 

Member  force  calculation 

- 

0.15 

6 

Analysis  uncertainty 

0.25 

17 

Total  uncertainty  in  ILF 

0.61 

Of  the  above  uncertainties,  those  due  to  natural  period,  damping,  member  force  calcu¬ 
lation,  and  analysis  uncertainty  can  be  considered  as  uncertainties  in  structural  analysis. 
Clearly,  the  uncertainty  in  estimating  the  natural  period  is  the  most  significant.  Indeed,  a 
small  error  in  this  quantity  leads  to  a  large  error  in  the  calculated  nominal  member  forces, 
and  stresses. 

Thus,  emphasis  should  be  given  in  improving  the  accuracy  of  methods  for  estimating 
the  natural  period  of  platforms.  This  could  be  achieved  by  improving  the  accuracy  of  pro¬ 
cedures  for  estimating  the  mass  of  a  platform,  by  modeling  the  behavior  of  the  foundation 
more  accurately,  and  by  accounting  for  the  joint  flexibility. 
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7.0  Fatigue 


Fatigue  is  an  important  consideration  in  structural  design.  For  many  structural  sys¬ 
tems,  fatigue  is  the  most  critical  failure  mode,  and  thus  safety  requirements  associated 
with  fatigue  reliability  dictate  design  decisions.  Fatigue  strength  can  be  described  by  a 
characteristic  S-N  curve  or  by  a  fracture  mechanics  fatigue  model. 

A  cumulative  damage  based  approach  for  fatigue  analysis  consists  of  the  following  steps: 

a)  modeling  the  loading  environment, 

b)  modeling  loads  exerted  by  the  environment  on  the  structure, 

c)  evaluation  of  nominal  loads, 

d)  evaluation  of  stresses  at  all  points  of  possible  crack  initiation  (stress  concentrations), 

e)  evaluation  of  cumulative  fatigue  damage  over  the  lifetime  of  the  structure. 

Uncertainties,  which  are  involved  in  steps  a)  to  c),  were  studied  in  the  previous  chapters. 
In  this  chapter,  we  quantify  errors  in  calculating  stress  concentration  factors.  We  also 
combine  the  errors  involved  in  all  the  steps  of  fatigue  analysis  and  quantify  uncertainties 
in  fatigue  damage,  for  the  case  that  a  cumulative  damage  approach  is  used.  Finally,  we 
investigate  the  relaf  1  r  oortance  of  the  uncertainties  in  each  of  the  steps  a)  to  e),  and 
also  of  the  random  and  modJiag  uncertainties. 

7.1  Uncertainties  in  stress  concentration  factor. 

Wirsching  (81)  reported  estimates  of  the  uncertainties  in  stress  concentration  factors 
for  tubular  joints  of  offshore  structures.  These  uncertainties  are  for  stress  concentration 
factors  obtained  from  parametric  equations,  such  as  those  by  Kuang,  Potvin  and  Leick 
(77).  According  to  Wirsching,  the  average  bias  in  the  stress  concentration  factor  is  in 
the  range  from  0.80  to  1.20,  and  the  COV  ranges  between  0.1  and  0.50.  The  bounds  for 
the  bias  and  the  COV  in  stress  concentration  factor  are  very  wide,  possibly  because  the 
parametric  equations  cover  a  large  number  of  geometries  and  loading  conditions. 

Uncertainties  in  the  stress  concentration  factor  axe  large  for  other  engineering  structures. 
For  example,  the  stress  concentration  factor  for  the  fatigue  analysis  of  a  liquid  propellant 
engine  was  assumed  to  follow  the  beta  distribution  (Sutharsana  et  al  90),  (Moore,  et  al 
90).  The  stress  concentration  factor  is  in  the  range  from  1.2  to  3.5,  and  its  COV  is  roughly 
0.15. 
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7.2  Uncertainties  in  cumulative  fatigue  damage. 

Studies  on  fatigue  reliability  of  marine  structures  assume  that  the  effect  of  random 
uncertainties  on  the  cumulative  damage  is  negligible.  Thus,  these  studies  account  only  for 
modeling  uncertainties  in  stress  evaluation  processes  (Wirsching  and  Chen  88).  Although 
the  effect  of  random  uncertainties  reduces  with  the  number  of  load  cycles  increasing,  to 
the  best  of  our  knowledge  no  study  has  proven  that  this  effect  is  negligible. 

The  objectives  of  the  exercise  presented  in  this  section  are  to  address  the  above  issue, 
estimate  uncertainties  in  the  cumulative  fatigue  damage  over  the  lifetime  of  platforms  and 
ships,  and  study  the  relative  importance  of  each  uncertainty. 


The  following  are  the  basic  assumptions: 

a)  Fatigue  life  can  be  estimated  by  using  the  S-N  curves.  The  integrated  Paris  law  gives 
the  same  form,  assuming  no  threshold  level  and  no  sequence  effects.  The  slope  of  these 
curves  is  constant  for  any  number  of  cycles,  N. 

b)  Miner’s  rule  can  be  used  to  estimate  fatigue  damage. 

c)  The  stress  amplitude  distribution  is  known. 

d)  The  mean  and  standard  deviation  of  the  cumulative  damage,  D,  can  be  estimated  by 
linearizing  the  expression  relating  D  with  ail  random  variables  around  the  mean  values 
of  these  variables.  This  is  a  crude  approximation  because  the  derivatives  of  the  damage 
with  respect  to  the  values  of  the  random  variables  are  not  constant.  Advanced  methods 
for  fast  probability  integration  are  more  accurate  in  this  case  (Madsen,  Krenk  and  Lind 
86).  However,  the  objective  of  this  study  is  to  identify  the  most  important  uncertainties 
and  obtain  only  rough  estimates  of  the  COV  of  D.  Moreover,  the  estimates  for  the 
bias  and  the  COV  of  the  random  variables  involved  in  damage  calculations  are  very 
crude.  Thus,  the  benefits  from  using  an  advanced  fast  probability  integration  method 
are  minimal. 


Under  the  above  assumptions,  fatigue  damage  can  be  calculated  by  the  following  equa¬ 
tion  (Wirsching  and  Chen  88): 


D  vnjjsr 

A 


(7.1) 


where, 


Bn  represents  the  modeling  error  in  the  stress  at  points  of  stress  concentration, 
m,  is  the  exponent  in  the  S-N  curves, 

Si ,  is  the  predicted  stress  amplitude  at  the  ith  load  application,  and 
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i4,  is  the  constant  at  the  right  hand  side  of  the  S-N  equations. 

The  summation  is  for  all  load  applications. 

The  above  equation  is  based  on  the  assumption  that  the  stress  is  a  narrowband  process. 
The  stress  might  be  wideband  which  affects  the  fatigue  damage  and  the  life  of  the  structure. 
We  may  account  for  the  above  effect  by  using  the  rainflow  correction  factor  (Wirsching  84, 
Wirsching  and  Chen  88). 

The  modeling  bias  Bn  is  given  by  the  following  equation: 


Bu  =  Bm  •  Bs  •  Bf  •  Bn  •  Bh  , 


(7.2) 


where, 

Bm  represents  uncertainties  in  the  geometry  due  to  manufacturing  imperfections, 

Bs  represents  uncertainties  in  seastate  description, 

Bf  represents  uncertainties  in  wave  load  predictions, 

Bn  is  the  bias  for  errors  in  static  and  dynamic  structural  analysis,  and 

Bjj  is  the  bias  for  uncertainties  in  stress  concentration  factors. 

By  using  a  first  order  Taylor  series  expansion  of  the  expression  for  D  about  the  mean 
values  of  all  random  variables,  we  obtain  the  mean  value  of  D, 


E(D)  = 


Em(Bn)  £  Em(Si) 
B{A) 


(7.3) 


Assuming  that  the  statistics  of  the  predicted  stress  are  the  same  for  all  load  cycles,  we 
have, 

E(D)  =  (7.4) 

where  N  is  the  number  of  cycles  over  the  lifetime  of  the  ship. 

The  coefficient  of  variation  of  fatigue  damage  D  in  (7.1)  is, 


COVD  =  (m2COVgn  +  COV\  +  COVgSm)1/2,  where  (7.5) 

COVb,,  is  the  COV  of  modeling  bias, 
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COV  a  is  the  COV  of  A,  and 
COVzsm  is  the  COV  of  the  sum  ES™. 


Note  that  subscript  i  has  been  dropped. 

The  first  term  in  the  expression  with  the  square  root,  on  the  right  side  of  (7.5),  represents 
the  effect  of  modeling  uncertainties.  The  second  term  is  associated  with  uncertainties  in 
S-N  curves  and  the  third  represents  the  effect  of  random  uncertainties. 

As  it  was  mentioned  earlier,  equations  (7.4)  and  (7.5)  are  approximate.  The  reasons  for 
using  them  have  been  mentioned  earlier  in  this  section. 

7.2.1  Relative  importance  of  random  uncertainties. 

All  studies  on  fatigue  reliability  of  marine  structures  neglect  the  effect  of  random  uncer¬ 
tainties  on  the  fatigue  damage  without  justifying  this  approximation.  Here,  we  compare 
the  effect  of  random  uncertainties  against  that  of  modeling  uncertainties  in  order  to  assess 
the  importance  of  random  uncertainties.  We  also  investigate  the  effect  of  the  correlation 
between  the  maxima  of  the  stress  process. 

We  considered  two  cases.  In  the  first  case  the  maxima  of  the  stress,  Si,  follow  the 
Rayleigh  distribution,  while  in  the  second  they  follow  the  Weibull  distribution.  We  assume 
that  the  correlation  coefficient  between  the  ith  and  the  kth  stress  maxima,  Si  and  Sk,  is, 

PSiSh  =  Ps4"s*+,»  where  pstsi+i  is  the  correlation  coefficient  between  two  subsequent 
peaks.  In  our  study  we  considered  different  values  for  PSiSi+i  ifl  the  range  from  0.  to  0.99. 
In  the  following  discussion  subscripts  will  be  dropped.  After  some  algebra,  the  following 
equation  was  derived  for  the  COV  of  random  uncertainties: 


COVZa 


COy»m({±f  -  g&py" 
N  >/» 


where  COVs  is  the  COV  of  a  local  maximum. 


(7.6) 


It  is  observed  from  (7.6),  that  the  COV  increases  with  the  correlation  coefficient  between 
subsequent  maxima  increasing.  Moreover,  the  COV  for  random  uncertainties  decreases, 
with  the  number  of  load  cycles,  N ,  increasing.  The  COV^am  is  almost  zero  for  large  values 
of  N  (say  107),  for  any  value  of  p  less  than  one. 

The  COV  for  random  uncertainties  is  presented  in  Table  7.1,  for  three  cases.  In  the  first 
case,  the  stress  amplitude  follows  the  Rayleigh  distribution,  while  in  the  latter  two  cases, 
it  follows  the  Weibull  distribution  (Ochi  90)  with  Weibull  exponent  c  equal  to  0.7  and  1.0, 
respectively. 
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Table  7.1  Coefficient  of  Variation  of  Cumulative  Fatigue  Damage 
Due  to  Random  Uncertainties  (N  =  106) 


Distribution  of  COV  of  Cumulative  Damage 

Stress  Amplitude 

Correlation  Coefficient  of  Subsequent  Peaks 


0.0 

0.5 

0.8 

0.9 

0.99 

Rayleigh 

0.0 

0.0 

0.01 

0.01 

0.03 

Weibull  (c  =  1.0) 

0.01 

0.01 

0.01 

0.02 

0.06 

Weibull  (c  =  0.7) 

0.01 

0.01 

0.02 

0.03 

0.10 

It  is  observed  that  the  effect  of  random  uncertainties  is  small.  Moreover,  a  similar 
calculation  for  N  =  108,  which  is  a  typical  number  of  load  applications  over  the  lifetime  of 
a  marine  structure,  showed  that  the  COV  due  to  random  uncertainties  is  practically  zero. 
Since  the  distributions  considered  for  the  stress  peaks  represent  real  life  situations,  we 
conclude  that  random  uncertainties  can  be  neglected  in  fatigue  reliability  analysis  without 
losing  any  accuracy.  This  is  true  even  for  the  case  for  which  the  adjacent  stress  maxima  are 
strongly  correlated.  Hence,  the  value  of  p  is  also  unimportant  provided  that  the  number 
of  load  cycles  is  large  (say  107). 

7.2.2  Relative  contribution  of  various  types  of  uncertainty  on  fatigue  damage. 

The  equations  of  the  previous  sections  allow  to  quantify  the  uncertainties  in  the  cumu¬ 
lative  fatigue  damage.  Eq.  (7.4)  can  be  used  to  estimate  the  average  bias  while  eq.  (7.5) 
is  for  the  COV.  We  used  the  above  equations  to  calculate  the  average  bias  and  COV  of 
fatigue  damage  for  typical  marine  structures.  We  also  studied  the  relative  contribution  of 
various  uncertainties  to  the  overall  uncertainty  in  fatigue  damage. 

The  data  on  various  uncertainties,  which  are  involved  in  fatigue  analysis,  are  presented 
in  Table  7.2.  The  COV’s  in  the  first  column  are  based  on  discussions  with  experts  from 
the  offshore  industry.  (Wirsching  81).  Bea’s  data  are  for  offshore  platforms  while  Kaplan’s 
data  are  for  ships.  Kaplan  (84)  and  Bea  (89)  provided  estimates  for  Bs  and  Bp  only. 
In  addition,  the  value  of  COVbn  for  ships  was  estimated  in  Chapter  6.  The  COV  of  the 
other  variables  were  assumed  to  be  equal  to  the  corresponding  values  specified  in  the  first 
column. 

Finally,  exponent  m  was  taken  to  be  4.38. 
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Table  7.2  Uncertainties  Involved  in  All  Steps  of  Fatigue  Analysis 


Type  of  Uncertainty  COV 

Bea  (89)  Kaplan  (84) 


Bm 

0.2 

0.2 

0.2 

Bs 

0.5 

0.58“ 

0.15 

Bp 

0.2 

0.196 

bn 

0.3 

0.3 

0.12 

bh 

0.3 

0.3 

0.3 

Bca 

1.0 

1.0 

1.0 

Notes: 

a)  This  COV  refers  to  the  cumulative  effect  of  environmental  and  load  evaluation  uncer¬ 
tainties,  i.e.  to  the  product  of  Bs  and  Bp. 

b)  This  COV  represents  modeling  uncertainties  in  the  combined  wave  and  slamming  bend¬ 
ing  moment.  The  estimate  was  based  on  COV’s  of  0.1  and  0.16  for  the  uncertainties  in 
wave  and  slamming  bending  moments,  respectively. 

c)  Ba  is  the  bias  for  the  constant  A  at  the  right  hand  side  of  the  S-N  equations. 


There  is  uncertainty  in  stress  calculations  due  to  the  fact  that  the  peaks  of  the  wave 
elevation  do  not  follow  the  Rayleigh  distribution.  This  uncertainty  introduces  conservative 
bias  in  the  stress  predictions. 

Table  7.3  presents  the  overall  uncertainty  in  cumulative  damage,  as  well  as,  the  relative 
contribution  of  all  types  of  uncertainties  to  the  cumulative  damage. 
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Table  7.3  Uncertainties  in  Cumulative  Fatigue  Damage  and  Relative 
Contribution  of  Each  Uncertainty  in  Table  7.2 


Source 

COVD 

COVM 

COVs 

COVF 

COVN 

COVH 

COVvs^COVa 

3.29 

0.07 

0.45 

0.07 

0.16 

0.16 

0.0 

0.09 

Bea  (89) 

3.42 

0.07 

0.55 

0.15 

0.15 

0.0 

0.08 

Kaplan  (84) 

2.21 

0.16 

0.09 

0.14 

0.06 

0.35 

0.0 

0.20 

The  following  conclusions  can  be  extracted  from  Table  7.3. 

a)  Uncertainty  in  cumulative  damage  is  very  large  for  both  ships  and  offshore  platforms. 
The  reason  is  that  fatigue  damage  is  extremely  sensitive  to  the  amplitude  of  the  applied 
stress.  In  other  words,  a  small  change  in  the  amplitude  results  to  a  large  change  in  the 
fatigue  damage  and  the  expected  fatigue  life. 

b)  As  indicated  by  our  example,  uncertainty  in  fatigue  damage  is  smaller  for  ships  than 
for  offshore  structures. 

c)  Uncertainty  in  describing  the  loading  environment  is  the  most  important  for  offshore 
platforms.  This  means  that  even  a  small  reduction  in  this  uncertainty  will  result  to  a 
large  reduction  in  the  overall  uncertainty  in  fatigue  damage. 

d)  For  the  case  of  ships,  uncertainty  in  the  stress  concentration  factor  is  the  most  important. 
The  next  important  uncertainty  is  that  in  A,  which  is  the  constant  in  the  right  hand 
side  of  the  expression  for  the  S-N  curves. 

e)  The  effect  of  random  uncertainties  is  negligible  because  these  uncertainties  are  averaged 
out  in  the  procedure  for  evaluating  fatigue  damage.  Moreover,  the  statistical  correlation 
between  consecutive  stress  peaks  is  unimportant  in  fatigue. 
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8.0  Conclusions 


The  following  are  the  main  conclusions  from  our  study. 

1.  Tn  both  ships  and  offshore  platforms,  it  is  important  to  distinguish  between  natural  (ran¬ 
dom)  and  modeling  uncertainties  and  between  their  effects  on  lifetime  extreme  loads  and 
load  effects.  Since  the  latter  uncertainties  are  systematic,  they  do  not  decrease  with  the 
length  of  the  exposure  period  increasing.  Therefore,  if  we  treat  modeling  uncertainties 
as  if  they  were  random,  we  may  grossly  underestimate  uncertainties  in  extreme  loads. 
In  our  opinion,  approach  1  in  Table  4.9  is  more  appropriate  than  approach  2. 

2.  We  believe  that  modeling  uncertainties  have  not  been  treated  correctly  (i.e.  according 
to  the  method  described  in  the  previous  paragraph)  in  the  new  Load  and  Resistance 
Factor  Design  code  of  the  American  Petroleum  Institute.  As  a  result,  the  coefficients  of 
variation  of  load  effects  have  been  grossly  underestimated. 

3.  An  effective  way  to  reduce  uncertainties  is  to  distinguish  between  different  classes  of 
ships  with  different  characteristics  and  operational  schedules,  and  quantify  uncertainties 
for  each  class,  separately. 

4.  Moreover,  we  may  consider  the  dependence  of  modeling  errors  on  some  parameters,  such 
as  the  significant  wave  height,  and  the  relative  heading  angle.  Then,  we  can  employ 
regression  to  determine  relations  between  the  error  and  the  above  parameters.  Such 
relations  can  be  used  to  improve  theoretical  predictions  of  loads  and  load  effects.  From 
work  dealing  with  this  approach,  it  was  demonstrated  that  we  can  reduce  uncertainties 
by  a  significant  amount. 

5.  There  is  significant  variability  in  Stillwater  load  effects  in  different  voyages  of  a  ship. 
There  is  also  significant  variation  between  similar  ships  and  between  different  ship  types. 

6.  Different  effects  lead  to  underestimation  of  overestimation  of  loads  or  load  effects  on 
ships.  For  example,  if  we  neglect  nonlinearities,  wave  directionality,  and  heavy  weather 
countermeasures,  then  we  will  underestimate  load  effects  in  ships.  On  the  other  hand, 
we  overestimate  load  effects,  if  we  do  not  account  for  errors  in  visual  observations  of 
wave  heights,  and  for  the  dependence  of  short-crestedness  on  sea  severity. 

7.  The  coefficient  of  variation  for  extreme  wave  moments  in  ships  is  roughly  0.20.  The  bias 
(ratio  of  actual  over  predicted  value)  is  greater  than  one,  which  means  that  theoretical 
estimates  of  wave  loads  are  lower  than  the  actual  values.  However,  the  magnitude  of 
this  exceedence  above  one  depends  upon  the  type  of  ship  as  well  as  consideration  of  the 
particular  type  of  bending  (ie.  sagging  or  hogging). 

8.  The  coefficient  of  variation  in  extreme  global  loads  (overturning  moments  and  base 
shear  forces)  on  offshore  platforms  ranges  between  0.60  to  0.90.  While  these  numbers 
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are  large,  and  their  precise  magnitude  may  be  questioned,  the  actual  values  are  expected 
to  be  considerably  larger  than  that  in  A.P.I.  code.  The  trend  indicated  by  these  numbers 
appears  to  be  proper. 

9.  A  large  uncertainty  is  involved  in  combining  slamming  and  wave  bending  moments  in 
ships.  This  problem  cannot  be  treated  properly  by  using  Turkstra’s  rule,  or  the  peak 
coincidence  approximation,  although  a  reasonable  value  appears  to  be  indicated  by  the 
square  root  of  sum  of  squares  approach. 

10.  No  reliable  validated  method  presently  exists  for  calculating  slamming  bending  mo¬ 
ments,  and  local  hydrodynamic  pressures  on  the  ship  hull  due  to  impact  with  waves. 
These  loads  are  very  significant  in  reliability  analysis  of  ships.  Therefore,  the  necessary 
resources  should  be  allocated  for  expanding  the  state  of  knowledge  on  slamming. 

11.  Linear  seakeeping  methods  cannot  estimate  hydrodynamic  pressures  on  the  ship  hull 
with  acceptable  accuracy.  This  is  particularly  true  in  the  vicinity  of  the  bow  and  the 
stern  of  the  ship  hull.  These  methods  are  more  effective  in  calculating  global  loads  and 
wave  bending  moments. 

12.  Nonlinear  effects  are  more  important  for  ships  with  small  block  coefficients,  such  as 
containerships,  than  for  ships  with  large  block  coefficients  such  as  tankers  and  bulk 
carriers. 

13.  The  probability  distribution  of  wave  peaks  is  not  exactly  Rayleigh,  and  subsequent 
peaks  are  statistically  dependent.  This  introduces  conservative  error  in  the  estimates 
for  extreme  wave  height.  The  bias  (as  defined  in  this  report)  due  to  this  effect  is  roughly 
0.93  and  the  coefficient  of  variation  is  0.08. 

14.  Random  uncertainties  are  unimportant  in  fatigue  reliability  analysis  of  both  ships  and 
offshore  platforms.  Moreover,  the  statistical  correlation  between  subsequent  wave  peaks 
is  also  unimportant. 

The  following  conclusions  refer  to  the  relative  importance  of  the  uncertainties  involved 
in  stress  analysis. 

15.  It  is  general  consensus  that,  in  offshore  structures,  the  uncertainty  in  describing  the 
loading  environment  is  the  most  important.  For  those  fixed  platforms,  for  which  dy¬ 
namic  effects  are  insignificant,  the  largest  part  of  this  uncertainty  is  due  to  errors  in 
estimating  the  long  term  maximum  wave  height.  The  uncertainty  in  estimating  the  iong 
term  severest  sea  condition  over  the  lifetime  of  the  platform  is  primarily  responsible  for 
the  above  errors. 

16.  Uncertainties  in  the  value  of  the  drag  coefficient  in  Morison  equation  are  also  important, 
especially  when  calculating  the  forces  on  small  cylindrical  members  of  offshore  platforms. 
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17.  Uncertainties  in  describing  the  loading  environment  (wave  height  and  period)  are  the 
most  important  in  fatigue  analysis  of  offshore  platforms.  Errors  in  stress  concentration 
factor  and  in  structural  analysis  follow  in  terms  of  relative  importance.  It  should  be 
mentioned  that  errors  in  structural  analysis  are  primarily  due  to  errors  in  estimating 
the  natural  period  of  a  platform. 

18.  The  following  are  the  most  important  uncertainties  in  stress  analysis  on  ships: 

-  uncertainties  in  slamming  bending  moments, 

-  spectral  shape  variability, 

-  errors  in  combining  wave  and  slamming  bending  moments, 

-  variability  in  Stillwater  load  effects, 

-  model  uncertainty  in  the  response  amplitude  operator  for  wave  bending  moments. 
Note  that  the  order  in  which  the  above  uncertainties  are  presented,  is  not  associated 
with  their  relative  importance. 

19.  Uncertainty  in  stress  concentration  factor  is  the  most  important  factor  in  fatigue  analysis 
of  ships. 
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9.0  Research  Plan  for  Future  Efforts  to  Reduce  Uncertainties,  in  Estimated 
Load  Effects  on  Marine  Structures 

On  the  basis  of  the  information  given  in  the  preceding  chapters  of  this  report,  future 
research  should  be  directed  toward  various  projects  that  would  have  benefit  in  reducing 
the  most  important  uncertainties  in  load  effects  on  marine  structures.  A  five  year  research 
program  plan  is  provided  in  the  present  section  which  is  based  on  the  relative  importance 
of  the  various  tasks,  while  also  considers  the  associated  factors  of  cost  and  risk. 

The  tasks  are  briefly  described,  together  with  their  relative  ranking,  the  rationale  for 
such  a  prioritizing,  and  the  basis  for  assigning  the  cost,  benefit,  and  risk  characterization 
to  each  task.  Initially,  a  listing  of  the  tasks  is  given  with  their  associated  characterization 
(and  related  comments)  in  a  summary  form,  which  is  then  followed  by  the  task  descriptions 
in  the  following  paragraphs.  The  topic  of  cost  is  considered  in  terms  of  the  required  man¬ 
hours  to  complete  the  particular  task,  with  the  following  associated  definitions: 

Costs 


Low  1,000  -  1,500  man-hours 

Medium  2,500  -  3,500  man-hours 

High  4,000  -  6,000  man-hours 

A  proposed  schedule  for  carrying  out  the  work  proposed  in  this  research  program  plan 
is  appended  at  the  end  of  this  section. 
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Table  9.1  List  of  Tasks  for  Reducing  Uncertainties 


Task 

Import¬ 

ance 

Cost 

Risk 

Comments 

1.  Impact  loads  on  ships 

+ 

o 

o 

(local  effects) 

2.  Slam  loads 

+ 

- 

- 

Should  be  initiated 

(hull  grinder  responses) 

3.  Extreme  impact  loads 

+ 

o 

o 

after  completing  1 
Should  be  initiated 

4.  Combination  of  loads-load 

+ 

0 

o 

after  partially 
completing  2 

effects 

5.  Nonlinearities  in 

+ 

- 

- 

Will  complement 

in  wave  loads 

6.  Nonlinear  systems  under 

+ 

o 

+ 

SR- 1304 
project 

random  excitation 

7.  Modeling  uncertainties 

+ 

+ 

+ 

8.  Breakiig  waves 

+ 

+ 

+ 

9.  Effects  of  routing,  heavy 

o 

o 

- 

weather  countermeasures 

10.  Torsional  loads  in  ships 

0 

0 

o 

11.  Ship-superstructure  interaction 

o 

o 

- 

Very  important  for 

12.  Hydrodynamic  loads  in 

+ 

0 

+ 

ships  with  large 
superstructures. 
(Naval,  passenger 
ships) 

offshore  platforms 

13.  Natural  period  of 

+ 

+ 

+ 

offshore  platforms 
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Task 


Risk 


Comments 


14.  Flexibility  of  joints 


Import-  Cost 
ance 

o 


15.  Development  of  designer  o  4- 

oriented  computer  code 
for  stress  analysis 


Important  in 
transverse  strength 
analysis  of 
bulk,  ore  carriers, 
and  tankers. 

Should  be  initiated 
after  completing 
most  of  the 
previous  tasks. 


Notation:  +  :  High 

o  :  Average,  medium 
-  :  Low 

Note:  Order  reflects  priority  ranking. 
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Task:  Impact  loads  on  ships  (local  effects) 

Importance:  + 

Cost:  o 
Risk:  o 
Priority:  1 

Objective:  Compile  state-of-the-art  knowledge  on  impact  hydrodynamic  loads,  and 
also  extend  existing  approaches  to  establish  reliable  methods  for  predicting  hydrodynamic 
pressures  on  ship  hulls  due  to  bottom  slamming,  and  bow  flare  impact. 

Impact  pressure  due  to  slamming  and  bow  flare  impact  is  important  in  design  of  high 
speed  vessels  such  as  containerships,  SWATHS,  patrol  vessels,  and  combatants.  Indeed, 
impact  loads  are  known  to  induce  significant  damage,  usually  in  the  form  of  excessive 
permanent  set,  to  the  front  panels.  Evaluation  of  impact  pressure  is  a  very  difficult  and 
complex  task,  which  might  require  to  take  into  account  the  interaction  between  the  defor¬ 
mation  of  the  hull  panels  and  the  surrounding  water.  For  this  reason,  our  ability  to  predict 
this  type  of  load  and  its  effects  is  limited.  A  comprehensive  review  of  various  approaches  to 
the  problem  of  analyzing  hydrodynamic  impact  loads  was  presented  in  ISSC  (88)  (commit¬ 
tee  1.2),  as  well  as  in  previous  reports  of  the  same  committee.  Most  theoretical  approaches 
to  the  problem  are  based  on  a  momentum  equation  which  relates  the  impact  force  to  the 
change  of  momentum  of  the  water  surrounding  the  hull  (Kaplan  87).  A  serious  problem 
associated  with  these  methods  is  due  to  flat  bottom  surfaces  which  sometimes  introduce 
some  subjective  considerations  into  the  analysis.  Kaplan  developed  an  approach  which  cir¬ 
cumvents  the  problem  by  use  of  a  quasi-three  dimensional  method  to  determine  the  total 
impact  force  on  large  flat  surfaces.  That  method  avoids  thj  nfinite  force  values  found  in 
two-dimensional  sections  with  flat  bottoms.  However,  no  associated  pressure  information 
was  developed  in  that  approach.  In  addition,  the  incorporation  of  the  treatment  of  flat 
portions  of  a  ship  hull  bottom  region  together  with  the  remaining  part  of  the  immersed 
ship  sections  could  also  have  to  be  accomplished. 

Empirical  models  which  presume  a  linear  relation  between  the  impact  pressure  and  the 
square  of  the  impact  velocity  have  been  extensively  used  by  Ochi  (73),  Mansour  (82),  and 
Ferro  (85).  It  should  be  noted  that  the  high  variability  in  the  value  of  the  coefficient  relating 
the  impact  pressure  to  the  square  of  impact  velocity,  (Mansour  82),  is  indicative  of  the 
uncertainty  associated  with  this  approach.  The  impact  pressure  has  also  been  measured 
from  tests  in  which  plates,  wedges,  or  cones  drop  in  the  water  (Takemoto  84),  or  from 
actual  ships  (Rask  86). 

In  this  task,  it  will  be  necessary  to  assess  the  state  of  knowledge  in  predicting  slamming 
induced  hull  pressures  and  many  existing  theoretical  and  empirical  approaches  to  produce 
a  reliable  method  for  predicting  these  loads.  The  work  should  also  investigate  the  effect 
of  structural  deformation  due  to  the  loads,  assess  its  significance,  and  possibly  the: class  of 
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ship  hulls  for  which  it  should  be  accounted  for. 

Clearly,  this  is  a  difficult  and  complex  task.  Therefore,  the  associated  cost  and  risk  might 
not  be  low.  However,  this  task  is  very  important  because  it  is  necessary  to  predict  impact 
loads  in  a  reliable  way,  in  order  to  design  safe  and  efficient  ships.  This  particular  task,  or 
a  somewhat  similar  type  program,  is  also  included  in  the  FY-90  and  FY-91  recommended 
Ship  Structure  Committee  research  programs  (pp.  36-37) 

Task:  Slam  loads  (hull  girder  responses) 

Importance:  + 

Cost:  - 
Risk:  - 
Priority:  2 

Objective:  Develop  a  methodology  for  evaluating  slam-induced  bending  moments,  and 
stresses  in  ships. 

As  it  was  demonstrated  in  this  report  for  high  sea  state  conditions,  slam-induced  stresses 
may  be  of  the  same  order  of  magnitude  as  those  due  to  longitudinal  wave  bending  moments. 
Furthermore,  as  has  been  demonstrated  by  Vulovich  et  al  (89),  bending  stresses  due  to 
slamming  are  significant  at  forward  locations  where  lower  stresses  are  normally  expected. 
As  such,  the  slamming-induced  stresses  play  an  important  role  in  design  of  certain  classes 
of  ships  (container  ships,  naval  vessels,  and  different  advanced  marine  vehicle  concepts). 

There  is  uncertainty  in  calculating  hull  girder  responses  primarily  due  to  the  difficulties 
in  predicting  the  slam  impact  local  loads  and  pressures  (as  described  in  Task  1  above). 
On  this  basis,  it  is  suggested  that  a  study  using  the  results  of  Task  1  be  undertaken  to 
improve  the  ability  to  calculate  slam-induced  hull  girder  stresses.  Hull  flexibility,  and  its 
influence  in  determining  vibratory  responses,  should  be  included  using  the  methods  of 
Kaplan  (72),  Vulovich  (89)  or  other  similar  approaches  that  represent  the  ship  structure 
as  a  prismatic  beam.  By  means  of  such  an  approach,  it  is  not  anticipated  that  there  will 
be  great  difficulties  in  response  determination  as  long  as  the  task  on  local  impact  forces 
(Task  1)  has  been  successfully  accomplished. 

Task:  Extreme  impact  loads. 

Importance:  + 

Cost:  o 
Risk:  o 
Priority:  3 

In  structural  reliability  evaluations,  it  is  not  sufficient  to  know  the  statistics  of  the  values 
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of  loads  or  load  effects  at  any  time  instant.  We  also  need  to  know  the  extreme  values  over 
certain  exposure  periods.  Therefore,  it  is  important  to  predict  the  extreme  values  of  impact 
loads  and  their  effects. 

The  problem  of  determining  the  extreme  value  statistics  of  impact  loads  is  considerably 
harder  than  the  corresponding  problem  for  low  frequency  wave  induced  loads,  because  of 
the  nonlinear  behavior  of  the  ship.  More  specifically,  even  if  the  wave  elevation  were  treated 
as  a  Gaussian  stationary  process,  the  whipping  stress  would  not  be  Gaussian.  Therefore, 
its  local  peaks  would  not  follow  the  Rayleigh  or  the  Rician  distribution. 

An  estimate  of  the  extreme  value  statistics  can  be  obtained  by  expressing  the  impact 
loads  or  load  effects  in  terms  of  a  second-order  Volterra  series  and  by  following  an  approach 
introduced  by  Naess  (85).  The  basic  idea  behind  this  approach,  which  uses  the  Kac- 
Siegert  (47)  technique,  is  to  expand  the  quadratic  impulse  response  function  in  terms  of 
its  eigenfunctions.  This  allows  to  represent  the  response  as  a  weighted  sum  of  squares 
of  known  stationary  processes,  and  to  evaluate  the  average  upcrossing  frequency  of  some 
known  threshold.  We  believe  that  this  approach  is  suitable  for  this  problem  because,  as  it 
has  been  observed  in  experiments,  the  hydrodynamic  impact  pressure  as  well  as  the  impact 
force  are  proportional  to  the  square  of  the  impact  velocity. 

For  the  reasons  explained  above,  we  believe  that  this  is  an  important  task.  Further¬ 
more,  it  requires  the  adaptation  and  extension  of  a  highly  sophisticated  approach,  which 
is  reasonably  well  understood.  On  this  basis  we  expect  that  the  cost  required  is  medium. 
There  is  also  some  risk  that  the  resulting  extreme  value  predictions  might  not  correlate 
well  with  measurements,  requiring  the  development  of  an  alternate  method. 

Task:  Combination  of  loads  and  load  effects. 

Importance:  + 

Cost:  o 

Risk:  o 

Priority:  4 

Objective:  Improve  our  ability  to  statistically  combine  loads  and  their  effects  in  ships, 
and  to  predict  their  extreme  values. 

In  many  cases,  two  or  more  loads  or  load  effects  act  simultaneously  on  the  same  struc¬ 
tural  element.  Therefore,  in  probabilistic  analysis  and  design,  we  must  superimpose  these 
loads  and  estimate  the  extreme  values  of  the  combined  load  process.  Examples  of  loads  or 
load  effects  that  must  be  combined  are  vertical  and  horizontal  bending  moments,  Stillwa¬ 
ter  and  wave  bending  moments,  slamming  and  wave  induced  bending  moments,  and  axial 
forces  and  lateral  wave  pressure  applied  to  hull  plates. 


66 


In  general,  the  local  maxima  of  the  loads  to  be  combined  do  not  occur  at  the  came  time 
mstant.  Furthermore,  the  loads  might  be  correlated  or  even  nonstationary.  This  is  the 
ease  lor  slamming  and  wave- induced  bending  moments.  Therefore,  the  load  combination 
problem  is  very  complex,  and  for  some  cases  it  has  not  been  solved  in  a  satisfactory  way. 

Some  simple  procedures  for  combining  loads  and  load  effects  are  Turkstra’s  rule  (l  urk 
i:t.ra  70),  the  square  root  of  sum  of  squares  law  (Mattu  80),  and  an  approach  which  assumes 
that  the  load  processes  coincide.  The  aforementioned  approaches  are  fairly  straightfor- 
vard  and  as  such  they  are  popular  to  practioners.  They  have  also  been  used  by  several 
-archers.  Turkstra’s  rule  has  been  used  by  Ferro  (85)  to  combine  slamming  and  wave 
..luuced  bending  moments  in  ship  hulls.  Kaplan  (84)  has  also  used  the  square  root  sum 
of  squares  law  to  solve  the  same  problem.  Unfortunately,  these  approaches  are  based  on 
some  very  strong  assumptions,  which  in  many  practical  situations  are  unrealistic.  Turk¬ 
stra’s  rule,  for  example,  is  accurate  only  for  the  case  that  the  processes  to  be  combined 
•re  statistically  independent.  Similar  considerations  will  apply  to  the  square  root  of  sum 
<  >f  squares  relation. 

Afore  advanced  sophisticated  methods  are  also  available.  These  methods  calculate:  the 
uncrossing  rate  of  a  safe  domain  in  the  space  of  random  load  and  strength  variables  by 
h  vector  defined  by  the  loads  to  be  combined.  This  can  be  used  to  estimate  the  extreme 
value  statistics  of  the  combined  process  Madsen  et  al.  (86).  Bounds  for  the  uperossing  rate, 
i oil  are  inexpensive  to  calculate  can  be  found  by  the  point  crossing  method  introduced 
by  Larrabee  and  Cornell  (84).  This  method  was  extended  by  Winterstein  and  Cornell  (84). 
However,  the  point  crossing  method,  in  the  form  described  in  the  above  references,  cannot 
account  for  the  correlation  between  loads. 

Wen  (77),  and  Wen  and  Pierce  (83)  introduced  the  load  coincidence  method.  This 
■'■tl'cd  accounts  for  the  correlation  between  the  random  processes  to  be  combined  by 
’eclating  and  using  the  average  coincidence  rate.  It  also  estimates  the  average  crossing 
o  of  a  safe  region.  One  of  the  disadvantages  of  the  load  coincidence  method  is  that 
.->  not  directly  applicable  to  the  case  of  nonstationary  loads.  Moreover,  tlm  au-rnge 
ncidence  rate,  which  is  required  by  the  method,  is  not  always  known  or  easy  to  calculate. 

A.s  it  was  demonstrated  in  this  report,  the  accuracy  of  the  approach  employed  fe.  load 
•  .•rnbination  affects  significantly  the  estimates  of  the  extreme  loads.  Clearly,  i  ne  ^sue 
of  load  combination  is  critical  in  probabilistic  analysis  and  design.  Therefore,  it  is  worth 
undertaking  a  study  whose  objective  is  to  improve  the  ability  to  predict  the  extreme  values 
of  two  or  more  stochastic  load  or  load  effect  processes  acting  simultaneously  on  the  same 
structural  component.  This  will  be  accomplished  by  selecting,  adapting,  and  cxtoudTg 
existing  approaches  for  load  combination.  The  following  are  a  few  cases  that  must  b, 
analyzed: 

:)  low  frequency  wave  and  slamming  induced  longitudinal  ship  hull  bending  stresses, 
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b)  wave  induced,  and  Stillwater  bending  moments  in  ships, 

c)  horizontal,  vertical  and  torsional  wave  induced  bending  moment  in  ships, 

d)  dead,  live,  and  wave  loads  in  offshore  platforms. 

This  task  should  improve  our  ability  to  establish  the  statistics  of  extreme  loads  and 
design  loads  for  code  calibration.  For  this  reason,  this  task  is  important.  We  do  not 
anticipate  any  problems  in  completing  this  task.  Therefore,  the  associated  cost  and  risk 
are  low. 


Task:  Nonlinearities  in  wave  loads. 

Importance:  + 

Risk:  o 
Cost:  o 
Priority:  5 

Objective:  Review  methods  for  nonlinear  analysis  of  wave  loads  and  select  and  recom¬ 
mend  the  most  practical  ones. 

It  is  recognized  that  linear  strip  theory  gives  only  a  rough  approximation  of  low  fre¬ 
quency  wave  bending  moments,  and  that  its  results  might  be  misleading  when  it  is  applied 
to  estimate  load  in  extreme  sea  states.  Indeed,  nonlinear  effects  are  important  in  high  sea 
states,  in  which  case  linear  strip  theory  underestimates  wave  loads  (Jensen  81,  Vulovitch 
89). 

The  proposed  project  will  synthesize  the  bibliography,  including  the  results  of  the  SR- 
1304  project  on  strategies  for  nonlinear  analysis  of  marine  structures,  (sponsored  by  SSC) 
and  identify  the  most  practical  approaches  for  prediction  of  nonlinear  loads. 

The  following  sources  of  nonlinearity  should  be  considered, 

-  hydrostatic  and  hydrodynamic  effects  due  to  flaring  of  ship  sections, 

-  viscous  effects, 

-  large  rolling  motion,  which  introduces  coupling  between  vertical  and  horizontal  forces 
and  their  effects  on  loads, 

-  steep  waves. 

Perturbation  and  nonlinear  simulation  methods,  based  on  strip  theory,  should  he  con- 
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sidered  in  this  project.  Presently,  three-dimensional  methods  cannot  be  easily  applied  for 
load  prediction  in  practical  design  problems.  (ISSC  Comm.  1.2,  1988).  However,  there  is 
an  existing  Accelerated  Research  Program  on  Nonlinear  Ship  Motion,  that  is  supported  by 
the  U.S.  Navy  Office  of  Naval  Research  (ONR),  which  is  aimed  at  developing  techniques  for 
calculating  nonlinear  hydrodynamic  loads.  This  5  year  program  (which  started  recently) 
also  includes  model  testing  and  validation  efforts,  with  the  results  expected  to  provide 
some  useful  applications  for  the  present  area  of  structural  load  prediction  (without  any 
extensive  efforts  by  the  Ship  Structure  Committee  programs). 

This  is  an  important  task  because  nonlinearities  may  play  a  dominant  role  in  deter¬ 
mination  of  design  loads.  We  believe  that  the  project  can  be  accomplished  at  low  cost 
and  the  risk  associated  with  this  endeavor  (from  the  perspective  of  the  Ship  Structure 
Committee  programs)  is  also  low,  based  upon  the  efforts  being  carried  out  separately  in 
the  ONR  program. 

Task:  Nonlinear  systems  under  random  excitation. 

Importance:  + 

Risk:  + 

Cost:  o 

Priority:  6 

Objective:  Review  and  recommend  methods  for  random  vibration  analysis  of  nonlinear 
ship  responses. 

In  probabilistic  analysis  of  dynamic  systems,  we  need  to  derive  the  statistics  of  the 
response,  and  its  extreme  values  from  the  known  statistics  of  the  driving  exitation.  For 
the  case  of  nonlinear  systems,  this  is  a  very  challenging  problem  because  the  response  is 
non-gaussian  even  if  the  driving  excitation  is  a  stationary  gaussian  stochastic  process.  Due 
to  the  nonlinear  character  of  the  ship  response  in  high  sea  states,  where  the  exciting  forces 
are  expected  to  be  nonlinear  also,  it  is  very  important  to  determine  the  statistics  of  the 
response. 

This  project,  which  will  be  initiated  after  completing  the  task  on  nonlinear  loads,  will 
review  and  recommend  methods  for  determining  the  statistics  of  ship  responses,  loads,  and 
load  effects.  It  should  also  review  the  results  from  a  pertinent  project  sponsored  by  the 
Ship  Structure  Committee  (Chakrabarti  87)  in  terms  of  their  utility. 

The  following  methods  should  be  studied: 

-  equivalent  linearization  techniques, 

-  Markov  based  approaches, 
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-  methods  based  on  Volterra  series  expansion  of  the  response, 

-  Monte  Carlo  time  domain  simulation. 

The  investigator  of  this  project  must  have  a  good  understanding  of  random  vibrations 
and  level  crossing  problems.  As  a  consequence,  the  associated  cost  might  not  be  low.  Some 
problems  may  also  be  encountered  if  the  analytical  results  do  not  correlate  well  with  those 
from  Monte  Carlo  simulation  or  measurements.  Furthermore,  the  large  computational  cost 
required  by  most  of  the  aforementioned  methods  might  be  a  problem. 


'Task:  Modeling  uncertainties. 

Importance:  + 

Cost:  + 

Risk:  + 

Priority:  7 

Objective:  Develop  a  philosphophically  consistent  and  practical  method  to  model  quan¬ 
tify,  and  account  for  modeling  uncertainties  in  reliability  analysis,  design,  and  code  cali¬ 
bration. 

It  is  important  to  distinguish  between  natural  and  modeling  uncertainties.  Indeed, 
various  studies  on  reliability  analysis  and  design  of  marine  structures  yield  erroneous  results 
because  they  treat  modeling  uncertainties  incorrectly  (Bea  89,  Olufsen  90).  Unfortunately, 
there  is  no  consensus  on  how  to  distinguish  modeling  uncertainties  from  natural  ones,  and 
how  to  treat  them  properly. 

We  should  study  modeling  uncertainties  closer  and  in  more  detail.  The  proposed  project 
should  address  the  following  issues: 

a)  should  we  take  modeling  uncertainties  into  account  in  reliability  analysis,  design,  and 
design  code  calibration? 

b)  what  is  the  best  way  to  account  for  modeling  uncertainties? 

c)  how  do  the  errors  in  various  steps  of  stress  analysis  procedures  propagate  when  esti¬ 
mating  Ufetime  extreme  loads  and  their  effects? 

d)  are  the  existing  databases  on  measured  responses  of  marine  structures  sufficient  to 
quantify  modeling  uncertainties? 

e)  what  experiments  should  we  conduct,  and  what  measurements  should  we  take,  in  order 
to  quantify  modeling  uncertainties? 
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Clearly,  it  is  important  to  answer  the  above  fundamental  questions  before  delivering  a 
reliability  based  design  code  to  practitioners.  Otherwise,  we  take  the  risk  of  developing  a 
code  which  prescribes  designs  with  highly  inconsistent  safety  levels  (Bea  89,  Olufsen  90). 
For  this  reason,  this  task  is  important.  On  the  other  hand,  this  is  obviously  very  difficult 
and  it  requires  significant  resources.  Also,  due  to  the  philoshophical  nature  of  the  questions 
related  with  this  task,  and  the  poor  understanding  of  the  nature  of  modeling  uncertainties, 
there  is  a  high  risk  that  this  task  may  not  be  completed  in  a  satisfactory  way. 


Task:  Breaking  waves 

Importance:  + 

Risk:  + 

Cost:  + 

Priority:  8 

Objective:  Study  forces  due  to  breaking  waves. 

Breaking  waves  can  cause  significant  impact  loads  in  both  ships  and  offshore  structures. 
Unfortunately,  our  ability  to  predict  these  loads  is  limited.  (ISSC  85,  Committee  1.2). 
Ochi  and  Tsai  (84)  published  results  for  the  case  of  vertical  walls  of  cylinders  in  breaking 
waves.  Other  related  work  of  this  nature  was  also  carried  out  by  Kjeldsen  and  Dean  (86), 
where  the  forces  measured  in  breaking  waves  were  much  larger  than  for  non  breaking  wave 
conditions. 

This  task  will  critically  review  the  literature  in  the  area  of  breaking  waves  and  it  will 
provide  methods  for  predicting  the  resulting  loads.  Data  for  structures  other  than  cylinders 
would  have  to  be  obtained  for  application  to  a  larger  range  of  marine  structures. 

This  is  an  important  task  because  of  the  loads  associated  with  breaking  waves  are  large 
and  they  might  dominate  in  design.  Due  to  the  complexity  of  the  nature  of  breaking 
waves,  and  the  poor  understanding  of  the  associated  phenomenon,  the  cost  is  likely  to  be 
high  and  significant  difficulties  are  expected  in  this  study.  Experimental  measurements  in 
controlled  laboratory  tests  would  be  a  main  source  of  data  for  configurations  not  considered 
in  previous  work,  which  is  another  reason  for  expecting  a  high  cost. 
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Task:  Effects  of  routing,  heavy  weather  countermeasures,  and  cargo  load  dis¬ 
tribution  in  design 

Importance:  o 

Risk:  - 

Cost:  o 

Priority:  9 

Objective:  Develop  a  method  for  predicting  extreme  loads  and  load  effects  in  ships, 
which  accounts  for  the  effects  of  routing,  heavy  weather  countermeasures  and  cargo  load 
distribution  in  design. 

The  statistics  of  the  lifetime  extreme  loads  and  their  effects  can  be  obtained  by  com¬ 
bining  the  short  term  statistics  at  all  environmental  and  operational  states  weighted  by 
the  corresponding  probability  of  occurrence  of  each  state.  In  principle,  this  is  a  simple 
operation.  However,  the  long  term  extreme  loads  are  affected  by  many  factors  such  as 
countermeasures  taken  under  extreme  weather  conditions,  weather  routing,  load  distri¬ 
bution,  and  navigation  errors.  Petrie  (86),  Guedes  Soares  (90),  and  Hutchinson  (86), 
investigated  the  effects  of  some  of  the  above  factors.  Their  results  demonstrated  that  it  is 
important  to  account  for  these  factors,  since  these  factors  can  increase  critical  load  effects 
such  as  wave  bending  moments  (Guedes  Soares  90). 

The  objective  of  the  proposed  project  is  to  review  existing  approaches  for  taking  into 
account  the  effects  of  the  aforementioned  factors  and  integrate  them  into  procedures  for 
predicting  long  term  extreme  values  of  loads  and  their  effects. 

This  objective  can  be  accomplished  by  using  Guedes  Soares  (90)  approach.  The  key 
idea  is  to  derive  the  conditional  probability  density  function  of  operational  parameters, 
such  as  relative  ship  heading,  conditioned  upon  environmental  parameters,  such  as  the 
significant  wave  height,  by  simulating  operational  decisions  by  Monte  Carlo  simulation. 

We  believe  that  this  task  can  be  accomplished  at  moderate  cost  and  we  do  not  anticipate 
any  problems  or  risks  in  completing  it. 


Task:  Torsional  loads  in  ships 

Importance:  o 
Risk:  o 
Cost:  o 
Priority:  10 

Objective:  Develop  methods  for  predicting  torsional  loads  and  their  effects.  Torsional 
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loads  are  important  for  design  of  ships  with  large  hatch  openings  such  as  containerships. 
The  uncertainty  involved  in  estimating  torsional  moments  is  considerably  higher  than 
that  in  longitudinal  bending  moments.  The  reason  is  that  the  dynamic  system,  whose 
input  is  the  wave  elevation  and  output  the  resulting  torsional  loads,  involves  considerable 
nonlinearities  associated  with  roll  motion. 

A  study  should  be  undertaken  in  order  to  reduce  uncertainties  in  estimating  torsional 
loads.  This  effort  should  allow  us  to  predict  torsional  loads  their  statistics  as  well  as  their 
statistical  correlation  with  the  vertical  and  horizontal  bending  moments.  It  is  important 
to  know  this  correlation  in  order  to  be  able  to  estimate  the  Von  Mises  stresses  at  critical 
locations  in  the  deck. 

The  approach  we  propose  to  improve  our  ability  to  predict  torsional  loads  will  involve 
along  the  following, 

a)  more  precise  determination  of  line  of  action  of  lateral  shear  forces, 

b)  inclusion  of  nonlinear  roll  motion  terms,  in  equivalent  linearized  form,  in  load  equations, 

c)  more  precise  determination  of  distribution  of  vertical  CG  of  distributed  mass  of  analyzed 
ship. 

The  above  actions  are  consistent  with  present  analysis  methods,  and  can  be  accom¬ 
plished  by  efforts  that  require  a  moderate  (average)  cost.  The  expected  risk  is  also  mod¬ 
erate,  in  view  of  prior  experience  where  torsional  moments  have  not  been  predicted  with 
sufficient  accuracy. 


Task:  Ship-superstructure  interaction 

Importance:  o 
Risk:  - 
Cost:  o 
Priority:  11 

Objective:  Improve  accuracy  of  structural  analysis  procedures  by  taking  into  account 
the  interaction  between  the  ship  hull  and  the  superstructure. 

In  analysis  of  ships  with  long  superstructures,  the  interaction  of  the  latter  with  the 
ship  hull  is  important.  Ordinarily,  analysis  methods  do  not  model  the  superstructure  and 
neglect  its  effects  because  it  is  believed  that  this  is  a  conservative  idealization.  However, 
this  is  not  true  for  all  cases.  Large  stresses  and  resulting  cracks  have  been  observed  at 
locations  in  the  hull,  which  are  near  the  superstructure.  These  have  been  attributed  to 
sudden  changes  in  the  geometry  and  the  flexibility  along  the  ship  hull.  These  sudden 
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changes  are  due  to  the  presence  of  the  superstructure. 

For  this  reason  we  propose  to  modify  and  extent  computer  codes  for  structural  analysis 
so  that  the  superstructure  is  analyzed  together  with  the  hull. 

This  should  be  a  relatively  simple  task  which  should  also  improve  the  effectiveness  of 
stress  analyses  procedures.  The  associated  cost  and  risk  are  deemed  to  be  low. 

Task:  Hydrodynamic  loads  in  offshore  platforms 

Importance:  + 

Risk:  + 

Cost:  o 

Priority:  12 

Objective:  Improve  accuracy  of  load  calculation  methods  in  offshore  platforms. 

One  of  the  conclusions  from  this  study  is  that  the  uncertainty  in  hydrodynamic  loads 
in  fixed  offshore  platforms  is  the  most  important  of  all  loads  acting  on  such  structures. 
Therefore,  research  efforts  should  be  directed  toward  reducing  this  type  of  uncertainty. 

For  fixed  jacket  platforms,  Morison’s  equation  is  usually  employed  to  find  loads.  We 
should  improve  our  ability  of  predicting  wave  particle  kinematics,  and  of  estimating  the 
values  of  the  coefficients  in  Morison’s  equation.  We  should  also  investigate  alternative 
ways  to  estimate  loads.  For  floating  structures  first  order  diffraction  approaches  should 
be  further  validated  via  comparisons  with  test  data.  Second  order  diffraction  approaches 
should  also  be  further  developed  and  their  results  compared  with  measurements  particu¬ 
larly  for  low  frequency  drift  forces  and  also  high  frequency  springing  effection  TLP  and 
other  moored  structures. 

The  proposed  project  will  address  the  above  issues.  This  project  is  very  important 
because  uncertainties  in  hydrodynamic  loads  in  offshore  structures  are  critical.  This  should 
be  a  long  term  effort  and  the  required  cost  could  be  considerable.  The  risk  is  also  high  due 
to  the  extreme  complexity  and  length  of  the  effort.  However,  even  a  small  improvement  in 
the  accuracy  of  the  load  calculation  procedures  should  lead  to  significantly  more  efficient 
designs  (ISSC  85,  Committee  V2).  In  view  of  extensive  efforts  by  the  offshore  engineering 
industry,  most  of  the  work  will  be  accomplished  by  these  private  sources  rather  than 
through  Ship  Structure  Committee  studies. 
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Task:  Natural  period  of  offshore  platforms 

Importance:  + 

Risk:  + 

Cost:  4- 

Priority:  13 

Objective:  Reduce  the  error  in  procedures  for  calculating  the  natural  period  of  deep 
water  fixed  offshore  platforms. 

The  uncertainty  in  the  natural  period  of  fixed  offshore  platforms  is  the  most  critical 
in  dynamic  analysis.  (Moses  83).  This  uncertainty  is  primarily  due  to  variations  in  the 
platform  weight  and  the  added  hydrodynamic  mass,  and  the  uncertainty  in  the  platform 
stiffness.  The  latter  is  due  to  the  flexibility  of  the  joints  and  errors  in  modelling  the 
foundation.  Therefore,  we  should  reduce  this  critical  uncertainty  by  focusing  on  the  above 
sources  of  variability. 

The  cost  required  to  complete  this  task  might  be  high  from  the  point  of  view  of  total 
effort  devoted  to  this  type  of  study.  The  reason  is  that  it  will  be  necessary  to  develop  a 
model  for  the  behavior  of  flexible  joints  and  the  behavior  of  the  foundation  which  might 
be  very  complex  and  nonlinear.  It  is  sometimes  also  difficult  to  estimate  accurately  the 
added  mass  of  the  platform,  depending  upon  its  specific  design.  However,  in  view  of  the 
efforts  of  the  offshore  industry,  such  work  would  be  accomplished  primarily  by  these  private 
sources  rather  than  via  Ship  Structure  Committee  studies.  The  risk  in  this  area  may  be 
small,  since  a  design  allowance  for  variability  in  amplification  ratios  is  usually  considered 
in  platform  design. 

Task:  Flexibility  of  joints  and  rigid  length  of  transverse  webs  in  ships 

Importance:  o 

Risk:  - 

Cost:  - 

Priority:  14 

Objective:  Extent  computer  codes  for  ship  structural  analysis  so  that  the  joint  flexibility 
and  the  effect  of  rigid  length  of  transverse  webs  is  accounted  for. 

Hughes  (83)  (pp.  302-315),  demonstrated  the  importance  of  taking  into  account  the 
joint  flexibility  and  the  rigid  behavior  of  portions  of  transverse  webs  in  ship  hulls.  He  also 
established  simple  finite  element  models  which  account  for  these  effects.  However,  these 
models  have  not  incorporated  in  any  finite  element  code  for  ship  structural  analysis.  As  a 
result,  existing  codes  which  treat  joints  as  rigid  might  and  frames  as  flexible  all  over  their 
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length  result  in  considerable  errors  in  transverse  stress  analysis. 

We  believe  that  the  accuracy  of  available  codes  (MAESTRO,  DAISY)  can  be  signif¬ 
icantly  improved  by  implementing  the  method  developed  by  Hughes.  Furthermore,  the 
increase  in  the  complexity  of  the  analysis,  and  the  computer  cost  should  be  low.  The 
improved  codes  should  be  particularly  effective  for  analyzing  commercial  ships  such  as 
tankers,  ore  carriers,  and  bulk  carriers,  because  transverse  strength  plays  an  important 
role  in  their  design.  The  cost  and  risk  associated  with  this  task  are  low  because  the 
necessary  finite  element  models  and  theoretical  derivations  have  been  already  completed. 


Task:  Development  of  designer  oriented  computer  code  for  stress  analysis 

Importance:  o 

Risk:  - 

Cost:  + 

Priority:  15 

Objective:  Combine  computer  codes  on  seakeeping,  hydrodynamic  load  evaluation, 
structural  analysis,  load  combination,  analysis  of  long  term  extreme  loads,  and  fatigue 
analysis,  to  produce  a  design  oriented  integrated  computer  package. 

The  completion  of  some  of  the  previous  tasks  will  result  in  significant  advancements  in 
the  state  of  knowledge  in  the  area  of  ship  strength.  We  should  translate  these  advancements 
into  improvement  in  the  effectiveness  and  the  accuracy  of  stress  analysis.  Furthermore, 
we  should  make  the  methods  developed  in  the  previous  tasks  available  to  practitioners. 
Therefore,  we  propose  to  initiate  a  project,  which  will  integrate  the  methods  and  the 
software  developed,  to  produce  an  integrated  computer  package  for  stress  analysis  in  ships. 
Existing  computer  codes  for  seakeeping,  such  as  SCORES  or  ABS/SHIPMOTION,  and 
for  structural  analysis,  such  as  MAESTRO  or  DAISY,  can  be  used.  The  final  product  will 
be  a  fully  documented,  user  friendly  computer  program  which  can  be  used  by  practioners 
in  ship  design. 

We  believe  that  this  is  an  important  task.  The  cost  associated  with  the  task  is  expected 
to  be  considerable  but  not  excessive.  We  do  not  expect  any  risks  associated  with  this 
project. 
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Schedule 

Proposed  5  year  Research  Program  -  Reduction  of  Uncertainties  in  Load  Effects 
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Appendix  A:  Calculation  of  Uncertainties  in  Lifetime  Maximum  Loads  or  Load 
Effects 


The  bias  B,  of  the  annual  maximum  load  has  a  lognormal  probability  distribution. 
Therefore, 

(£nb  -  A)2 
-e  2a2 


1 


/b(6)  -  ffV® 


where 

<r  =  (fa(l  +  COVg))1'2, 

A  =  £n(EB)  -y, 

and  EB  and  COVb  are  the  mean  and  the  coefficient  of  variation  of  B,  respectively.  The 
lognormal  distribution  belongs  to  the  exponential  class  of  probability  distributions  because 
it  satisfies  von  Mises’  condition  (Ochi  90).  Therefore,  the  maximum  value  of  B,  over  an  N 
year  period  B^N\  follows  the  Type  I  asymptotic  extreme  value  probability  distribution, 


FB{n)  (b)  =  exp  (-e~a"(b  “  S**)) 


where,  Bn  is  the  most  probable  maximum  over  the  N  year  period,  and 


a.N  —  N  •  /b(Bn). 


The  most  probable  maximum,  Bjv,  satisfies  the  following  equation, 

P(B  >  Bn)  =  i, 


which  is  equivalent  to, 

P(CnB  >  £tiBn)  = 

Therefore,  since  £nB  is  normally  distributed  with  mean  A  and  standard  deviation  a, 

Bn  =  exp  ($~1(— )g  +  A) 

where  $(•)  denotes  the  probability  distribution  function  of  a  standard  Gaussian  random 
variable.  Finally,  the  coefficient  of  variation  of  B W  is 

u  w/V® 

VB(N)  =  - 5 - ; - 

QnBn  +  7 

where  a  is  the  Euler’s  constant  (0.577). 
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Appendix  B:  Combination  of  Slamming  and  Wave 
Induced  Stresses:  A  Simulation  Study 


B.l  Introduction 

The  objective  is  to  investigate  the  uncertainty  of  state-of-the-art  methods  for  superim¬ 
posing  slamming,  and  wave  induced  stresses.  For  this  purpose,  we  conducted  a  Monte- 
Carlo  simulation  study  to  compare  the  maximum  combined  slamming  and  wave  induced 
bending  moment  against  that  estimated  by  three  approximate  methods  for  load  combina¬ 
tion.  The  first  method  is  Turkstra’s  rule  and  it  has  been  utilized  by  Ferro  and  Mansour 
(85).  The  second  method  assumes  that  the  maxima  of  the  two  loads  coincide.  This  method 
will  be  called  peak  coincidence  approximation.  The  last  method,  which  is  called  square 
root  of  sum  of  squares  (SRSS),  has  been  suggested  by  Mattu  (80).  The  SRSS  method  has 
also  been  employed  by  Kaplan  (84)  to  combine  slamming  and  wave  bending  moments. 

Other  generic  methods  for  combining  loads  are  the  load  coincidence  (LC)  method,  which 
was  introduced  by  Wen  (77,  82),  and  the  point  crossing  method  proposed  by  Winterstein 
and  Cornell  (84).  LC  method  is  more  advanced  than  the  first  three  methods,  which  were 
mentioned  in  the  previous  paragraph,  and  it  accounts  for  the  statistical  dependency  of 
the  load  processes  to  be  combined.  However,  it  requires  some  information  on  the  load 
processes  which  is  difficult  to  be  estimated  for  the  case  of  slamming  and  wave  bending 
combination.  For  example,  it  is  very  difficult  to  calculate  the  mean  rate  of  coincidence 
of  the  aforementioned  loads,  or  the  conditional  probability  distribution  of  the  maxima  of 
the  combined  load  process,  given  that  the  maxima  coincide.  Therefore,  LC  method  is  not 
readily  applicable  to  the  problem  of  combining  bending  moments.  Moreover,  to  the  best 
of  our  knowledge,  it  has  never  been  applied  to  the  problem  of  slamming  and  wave  bending 
combination.  Therefore,  we  did  not  consider  this  method  in  this  study.  On  the  other  hand, 
Winterstein’s  method  assumes  that  the  load  processes  are  statistically  independent,  which 
is  not  a  realistic  assumption.  We  also  did  not  consider  this  method  in  our  study. 

The  procedure,  which  was  used  in  this  study,  consists  of  the  following  steps: 

a)  Select  a  representative  ship. 

b)  Select  a  ship  speed  and  a  significant  wave  height. 

c)  Simulate  a  time  history  of  the  wave  elevation,  by  using  a  random  number  generator. 

d)  Evaluate  the  resulting  wave  bending  time  history  using  the  SCORES  seakeeping  com¬ 
puter  code.  Find  the  extreme  value. 

e)  For  the  time  history  which  was  generated  in  step  c,  calculate  the  time  history  of  the 
resulting  slamming  bending  moment.  Find  the  extreme  value. 
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f)  Superimpose  the  wave  and  slamming  bending  moments,  which  were  generated  in  tasks 
d  and  e,  to  obtain  the  combined  bending  moment. 

g)  Find  the  extreme  value  of  the  combined  bending  moment.  Consider  this  value  as  the 
actual  maximum  combined  bending  moment. 

h)  Use  Turkstra’s  rule,  the  peak  coincidence  approximation,  and  the  S.R.S.S.  approach  to 
estimate  the  extreme  combined  bending  moment  from  the  time  histories  of  the  slamming 
and  wave  bending  moment,  which  were  generated  in  tasks  d  and  e. 

i)  Compare  the  estimated  combined  bending  moment  against  the  actual  bending  moment 
(task  h).  Seeps  c)  -  i)  are  repeated  in  order  to  obtain  a  sufficiently  large  number  of 
samples. 

The  data  used  in  this  study  are  presented  in  Table  B.l.  The  computer  program  SCORES 
was  employed  to  obtain  information  which  was  necessary  for  the  simulation,  such  as  two 
dimensional  properties,  and  response  amplitudes  and  phases  for  bending  moment,  relative 
motion,  and  velocity  as  function  of  frequency.  Random  time  histories  for  the  wave  elevation 
were  generated  according  to  the  procedure  suggested  by  Kaplan  (72,  87).  The  method 
described  by  Kaplan  (72,  86)  was  employed  to  calculate  the  slamming  bending  moment 
in  task  e.  Finally,  we  assumed  that  the  hull  is  rigid  in  calculating  the  slamming  bending 
moment  as  the  response  to  impact  forces.  This  assumption  might  lead  to  underestimation 
of  the  slamming  bending  moment  because  it  neglects  the  effects  of  the  hull  flexibility  and 
vibratory  responses.  However,  it  is  acceptable  in  the  context  of  our  study,  because  the 
objective  is  to  assess  the  error  in  combining  slamming  and  wave  bending  moments  and  not 
to  find  the  exact  values  of  these  quantities. 

Table  B.  1.  Data  for  Simulation  Study 


Ship  Length  (m): 
Displacement  (tons): 
Beam  (m): 

Draft  (m): 

Block  Coefficient: 
Speed  (knots): 


193 

36100 

27.6 

11.0 

0.6 

15,  18,  20,  25 


Type  of  Wave  Spectrum:  Pierson  Moskowitz  (1964) 
Significant  Wave  Height  (m):  6.1,  9.6,  13.1 
Length  of  Simulated  Time  Period  (min):  10 
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B.2.  Results  and  discussion 

The  estimates  of  the  maximum  combined  bending  moment  over  a  time  interval  of  10  min 
are  presented  in  Tables  B2-B7.  Results  from  Monte-Carlo  simulation,  Turkstra’s  rule,  peak 
coincidence  approximation,  and  the  SRSS  law  are  presented  in  the  above  tables.  These 
results  are  for  ship  speeds  from  15  to  25  knots  and  for  significant  wave  heights  in  the  range 
from  6  to  13  m.  The  bias  and  scatter  of  the  maximum  combined  bending  moment  are 
also  presented.  These  estimates  have  been  derived  from  200  simulated  time  histories  of 
the  combined  bending  moment. 

The  maximum  bending  moment  for  ship  speed  equal  to  15  knots  is  plotted  as  func¬ 
tion  of  the  significant  wave  height  in  Figure  Bl.  This  figure  compares  the  maximum 
simulated  bending  moment  against  that  estimated  by  Turkstra’s  rule,  peak  coincidence 
approximation,  and  SRSS  law.  The  shaded  regions  correspond  to  the  range  enclosed  by, 
a)  the  expected  maximum  combined  bending  moment,  and  b)  the  same  quantity  plus  one 
standard  deviation. 

The  characteristic  design  value  of  the  combined  bending  moment  at  a  risk  level  a,  is 
defined  as  the  threshold  which  corresponds  to  a  probability  of  at  least  one  upcrossing 
equal  to  a.  The  characteristic  design  values  for  probability  levels  0.05  and  0.01,  which 
were  obtained  by  simulation  and  by  the  two  approximate  methods  are  presented  in  Figure 
B2.  The  results  in  this  figure  correspond  to  a  ship  speed  of  15  knots  and  significant  wave 
height  ranging  from  6.14  m  to  13  m. 

It  is  observed  that  the  peak  coincidence  approximation  is  conservative  and  that  its  esti¬ 
mates  are  consistently  higher  than  the  values  obtained  by  simulation.  The  bias  of  the  peak 
coincidence  approximation  is  roughly  0.85  which  means  that  the  actual  combined  bending 
moment  is  15%  lower  than  the  estimated  value.  This  trend  of  the  peak  coincidence  ap¬ 
proximation  to  overestimate  loading  can  be  explained  as  follows.  The  maximum  slamming 
and  wave  bending  moment  do  not  occur  simultaneously.  This  can  be  observed  from  Figure 
B3  which  depicts  a  simulated  time  history  of  the  combined  bending  moment.  Apparently, 
the  maximum  of  the  combined  bending  moment  highly  depends  on  the  phase  lag,  S}  which 
corresponds  to  the  elapsed  time  between  the  occurence  of  a  slamming  and  a  wave  bending 
moment  peak.  The  most  hazardous  case  is  when  6  is  zero  which  means  that  the  peaks 
of  the  two  loadings  coincide.  It  is  known,  that  a  slam  almost  always  occurs  when  the 
bending  moment  changes  from  hogging  to  sagging.  Moreover,  phase  lag,  6 ,  is  not  likely  to 
be  zero  (Ochi  73).  Consequently,  the  peaks  of  the  two  loadings  are  not  likely  to  coincide. 
Therefore,  the  peak  coincidence  approximation,  which  is  based  on  the  aforementioned  un¬ 
realistic  assumption,  yields  conservative  results.  SRSS  law  ignores  the  trend  of  6  to  take 
values  within  a  narrow  range  and  therefore  its  estimates  are  also  in  error.  In  our  study, 
the  error  was  lower  than  that  from  the  peak  coincidence  approximation.  Furthermore,  the 
results  from  SRSS  law  did  not  show  any  trend  of  being  conservative  of  unconservative. 
This  observation  agrees  with  the  conclusions  of  Wen  (82). 
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In  general,  Turkstra’s  rule  is  not  conservative.  However,  its  results  were  conservative 
for  high  risk  levels.  The  large  scatter  of  the  estimates,  which  were  found  roughly  equal  to 
20-25%,  explains  this  observation.  The  histograms  of  the  maximum  combined  loading  from 
simulation  and  from  Turkstra’s  rule  (Figure  B4)  show  that,  although  the  average  value  of 
the  latter  is  lower  than  that  of  the  former,  the  estimates  from  Turkstra’s  rule  are  more 
likely  to  exceed  a  high  threshold  because,  the  corresponding  histogram  is  more  dispersed. 
This  relatively  high  scatter  of  Turkstra’s  rule  estimates  can  be  explained  by  the  equation 
defining  this  rule  for  the  case  of  two  load  processes: 

maxr{hr(0}  —  max^maxT{bw(t)}  -I-  b3(t),bw(t)  +  maxr{&a(t)}|, 

where  br(t),ba(t),  and  bw(t)  represent  the  combined  slamming,  and  wave  induced  bending 
moment,  and  maxr{-}  denotes  the  maximum  value  of  the  bracketed  process  over  a  time  pe¬ 
riod  of  length  T.  In  general,  the  coefficients  of  variation  of  random  variables  b3(t)  and  bw(t) 
are  considerably  higher  than  the  coefficients  of  variation  of  ma,XT{bw(t)} ,  maxT{bs{t)},  and 
maxT{br(t)}-  In  particular,  if  the  length  of  time  interval  T  is  large,  then  the  coefficients 
of  variation  of  the  maxima  are  very  low  compared  to  that  of  ba(t)  and  bw(t).  As  a  result, 
the  coefficient  of  variation  of  the  random  function  on  the  right  hand  side  of  Turkstra’s 
rule  equation  is  larger  than  that  of  the  actual  maximum  of  the  combined  stress  process, 
maxT{br(t)}. 

Finally,  the  effect  of  the  number  of  simulated  time  histories  on  the  estimated  bias  is 
shown  in  Figure  B5.  It  is  observed  that  the  number  simulated  time  histories  considered  in 
this  study  is  sufficient  for  estimating  the  above  statistics  with  good  accuracy. 

B.3.  Conclusions 

The  following  are  the  main  conclusions  from  our  study: 

a)  The  results  of  simplified  approaches,  which  are  used  to  combine  the  slamming  and 
wave  induced  bending  moments,  may  lead  to  serious  errors  in  the  estimated  design  loads. 
Note  that  the  bias  introduced  by  these  approaches  is  systematic,  and  as  such,  it  will  not 
decrease  even  if  the  length  of  the  time  interval,  which  is  considered,  becomes  very  high 
(say  20  years). 

b)  The  peak  coincidence  approximation  is  conservative  with  a  bias  in  the  range  from  0.60 
to  0.80. 

c)  The  results  from  SRSS  law  and  Turkstra’s  rule  do  not  consistently  lie  on  the  same 
side,  conservative  or  nonconservative,  of  the  actual  load  value.  Thus,  we  cannot  bracket 
the  expected  maximum  combmed  bending  moment  using  one  of  these  approaches  together 
with  the  peak  coincidence  approximation. 

d)  The  design  load  estimates  from  Turkstra’s  rule  have  significantly  larger  variability 
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Figure  Bl:  Maximum  Wave  Bending  Moment  vs.  Significant  Wave 
Height 

Note:  For  each  set  of  curves  lower  one  corresponds  to 

expected  value,  upper  one  corresponds  to  expected 
value  plus  one  standard  deviation. 
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Figure  B2:  Design  Bending  Moment  vs.  Significant  Wave  Height 
Note:  For  each  set  of  curves,  the  lower  one  corresponds  5% 
risk  level,  and  the  upper  one  to  1Z  risk  level. 
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Figure  BA:  Maximum  Combined  Bending  Moment  Histograms 

(The  histogram  traced  by  the  thick  line  corresponds  to 
simulation  results,  and  histogram  traced  by  the  thin  line 
corresponds  to  Turkstra’s  rule) 
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Figure  B.  5  Estimated  Bias  vs.  Number  of  Simulated  Time  Histories 


than  the  actual  load.  This  means  that  this  approximate  rule  introduces  both  bias  and 
variability,  which  is  added  to  the  statistical  variabil 


Table  B2.  Statistics  of  Maximum  Value 
of  Combined  Bending  Moment  for  Significant 
Wave  Height:  6.14m,  Ship  Speed:  15  knots 


Peak 

Method/  Turkstra’s  Coincidence 

Quantity _ Simulation _ Rule _ Rule _ SRSS 


Mean  (tn.m) 

0.8  xlO5 

Bias 

1.0 

Scatter 

0.09 

Char.  Design  Value 
(a  =  0.1)  (tn-m) 

0.89xl05 

Char.  Design  Value 
(a  =  0.05)  (tn-m) 

0.92xl05 

Char.  Design  Value 
(a  =  0.01)  (tn-m) 

0.96  x10s 

0.71x10s 

0.11x10® 

0.8x10 

1.1 

0.71 

1.0 

0.19 

0.15 

- 

0.86  x10s 

0.13x10® 

- 

0.95  x10s 

0.14x10® 

- 

0.13  x10s 

0.16x10® 

Average  Maximum  Wave  Bending  Moment  (tn.m)  =  0.648 xlO5 
Average  Maximum  Slamming  Bending  Moment  (tn.m)  =  0.472 xlO5 
Note:  All  bending  moments  correspond  to  sagging 
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Table  B3.  Statistics  of  Maximum  Value 
of  Combined  Bending  Moment  for  Significant 
Wave  Height:  9.6(m),  Ship  Speed:  15  knots 


Peak 


Method/ 

Quantity 

Simulation 

Turkstras 

Rule 

Coincidence 

Rule 

SRSS 

Mean  (tn.m) 

0.9xl05 

0.85x10s 

0.14x10® 

0.996xl05 

Bias 

1.0 

1.1 

0.64 

0.9 

Scatter 

0.143 

0.26 

0.162 

- 

Char.  Design  Value 
(a  =  0.1)  (tn-m) 

0.10  xlO6 

0.11x10® 

0.17x10® 

- 

Char.  Design  Value 
(a  =  0.05)  (tn-m) 

0.12  x10s 

0.13xl06 

0.18x10® 

- 

Char.  Design  Value 
(a  =  0.01)  (tn-m) 

0.13  xlO6 

0.17x10® 

0.21x10® 

- 

Average  Maximum  Wave  Bending  Moment  (tn.m)  =  0.72  xlO5 
Average  Maximum  Slamming  Bending  Moment  (tn.m)  =  0.695 xlO5 
Note:  All  bending  moments  correspond  to  sagging 
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Table  B4.  Statistics  of  Maximum  Value 
of  Combined  Bending  Mc.neni  :o.-  Significant 
Wave  Height:  13.07m,  Ship  Speed.  i  knots 


Method/ 

Quantity 

Simulation 

Turks  ras 
Rule 

Peak 

Coincidence 

Rule 

SRSS 

Mean  (tn.m) 

0.95  x10s 

0.91xlC5 

0.15x10® 

0.11x10® 

Bias 

1.0 

1.04 

0.62 

0.88 

Scatter 

0.15 

0.2u 

0.16 

- 

Char.  Design  Value 
(a  =  0.1)  (tn-m) 

0.11x10® 

o.i2xn® 

0.18x10® 

- 

Char.  Design  Value 
(a  =  0.05)  (tn-m) 

0.12x10® 

O.lSx.V3 

0.20x10® 

- 

Char.  Design  Value 
(a  —  0.01)  (tn-m) 

0.14x10® 

0.18x10° 

0.22  xl"3 

- 

Average  Maximum  Wave  Bending  Moment  (tn  n)  —  0.7a  xlO5 
Average  Maximum  Slamming  Bending  Moment  (tn.m)  =  3  77x  105 
Note:  All  bending  moments  correspond  to  sagging 
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Table  B5.  Statistics  of  Maximum  Value 
of  Combined  Bending  Moment  for  Significant 
Wave  Height:  6.14m,  Ship  Speed:  18  knots 


Peak 


Method/ 

Quantity 

Simulation 

Turkstras 

Rule 

Coincidence 

Rule 

SRSS 

Mean  (tn.m) 

0.1x10s 

0.83x10® 

0.14x10® 

0.99x10® 

Bias 

1.0 

1.2 

0.72 

1.0 

Scatter 

0.15 

0.21 

0.17 

- 

Char.  Design  Value 
(o  =  0.1)  (tn-m) 

0.12xl06 

0.11x10® 

0.17x10® 

- 

Char.  Design  Value 
(a  =  0.05)  (tn-m) 

0.14x10® 

0.12x10® 

0.19x10® 

- 

Char.  Design  Value 
(a  =  0.01)  (tn-m) 

0.15x10® 

0.15x10® 

0.2x10® 

- 

Average  Maximum  Wave  Bending  Moment  (tn.m)  =  0.7xl05 
Average  Maximum  Slamming  Bending  Moment  (tn.m)  =  0.71  xlO5 
Note:  All  bending  moments  correspond  to  sagging 
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Table  B6.  Statistics  of  Maximum  Value 
of  Combined  Bending  Moment  for  Significant 
Wave  Height:  6.14m,  Ship  Speed:  20  knots 


Peak 


Method/ 

Quantity 

Simulation 

Turkstras 

Rule 

Coincidence 

Rule 

SRSS 

Mean  (tn.m) 

0.12x10® 

0.98  x10s 

0.16x10® 

0.11x10® 

Bias 

1.0 

1.3 

0.76 

1.09 

Scatter 

0.22 

0.30 

0.19 

- 

Char.  Design  Value 
(a  =  0.1)  (tn-m) 

0.16x10® 

0.13x10® 

0.2x10® 

- 

Char.  Design  Value 
(a  =  0.05)  (tn-m) 

0.18x10® 

0.15x10® 

0.23x10® 

- 

Char.  Design  Value 
(a  =  0.01)  (tn-m) 

0.21x10® 

0.20x10® 

0.25x10® 

- 

Average  Maximum  Wave  Bending  Moment  (tn.m)  =  0.763 xlO5 
Average  Maximum  Slamming  Bending  Moment  (tn.m)  =  0.848 xlO5 
Note:  All  bending  moments  correspond  to  sagging 
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Table  B7.  Statistics  of  Maximum  Value 
of  Combined  Bending  Moment  for  Significant 
Wave  Height:  6.14m,  Ship  Speed:  25  knots 


Peak 


Method/ 

Quantity 

Simulation 

Turkstras 

Rule 

Coincidence 

Rule 

SRSS 

Mean  (tn.m) 

0.18x10® 

0.13x10® 

0.21x10® 

0.15x10® 

Bias 

1.0 

1.4 

0.83 

1.20 

Scatter 

0.21 

0.26 

0.16 

- 

Char.  Design  Value 
(a  =  0.1)  (tn-m) 

0.23x10® 

0.17x10® 

0.25x10® 

- 

Char.  Design  Value 
(a  =  0.05)  (tn-m) 

0.25x10® 

0.18x10® 

0.27x10® 

- 

Char.  Design  Value 
(a  =  0.01)  (tn-m) 

0.28x10® 

0.22x10® 

0.30x10® 

- 

Average  Maximum  Wave  Bending  Moment  (tn.m)  =  0.94 xlO5 
Average  Maximum  Slamming  Bending  Moment  (tn.m)  =  0.117xl06 
Note:  All  bending  moments  correspond  to  sagging 
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